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In a recent paper new upper and lower limits were given, in the context of the Schrodinger or 
Klein-Gordon equations, for the number No of S-wave bound states possessed by a monotonically 
nondecreasing central potential vanishing at infinity. In this paper these results are extended to the 
number Ne of bound states for the l-ih partial wave, and results are also obtained for potentials 
that are not monotonic and even somewhere positive. New results are also obtained for the case 
treated previously, including the remarkably neat lower limit Ne > {{[a/(2£ + 1) + 1] /2}} with 

a — (2/V) max r |V(r)| 1//2 (valid in the Schrodinger case, for a class of potentials that includes 

0<r<oc L J 

the monotonically nondecreasing ones), entailing the following lower limit for the total number N 
of bound states possessed by a monotonically nondecreasing central potential vanishing at infinity: 
N > {{(er + 1) /2}} {{(a + 3) /2}} /2 (here the double braces denote of course the integer part). 



I. INTRODUCTION AND MAIN RESULTS 

In a previous paper Q new upper and lower limits were provided for the number Nq of S-wave bound 
states possessed, in the framework of the Schrodinger or Klein-Gordon equations, by a central potential 
V(r) vanishing at infinity and having the property to yield a nowhere repulsive force, so that, for all 
(nonnegative) values of the radius r, 

V'(r) > 0, (1) 

hence 

-V(r) = \V(r)\. (2) 

In and always below, appended primes signify of course differentiation with respect to the radius r. 
The main purpose of the present paper is to extend the results of Q to higher partial waves, namely 
to provide new upper and lower limits for the number Ng of ^-wave bound states possessed by a central 
potential V(r). Here and always below i is of course the angular momentum quantum number (a non- 
negative integer). For simplicity we restrict attention here to the Schrodinger case, since the extension 
of the results to the Klein-Gordon case is essentially trivial, see Q . As in 0| we assume the potential to 
be finite for < r < oo, to vanish at infinity faster than the inverse square of r, 

lim \r 2+e V(r)] = 0, (3) 

r — >oo 

and, unless otherwise specified, not to diverge at the origin faster than the inverse square of r, 

lim [r 2 - £ V(r)] = 0. (4) 
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Here of course e denotes some positive quantity, e > 0. Moreover, in the following the "monotonicity" 
property Q is generally replaced by the less stringent condition 



or equivalently 



V(r) r 



[V{r)r- U ]' > 0, 



(5a) 



(5b) 



which is of course automatically satisfied by monotonic potentials vanishing at infinity, see and @ ; and 
we also obtain results for potentials that do not necessarily satisfy for all values of r the "monotonicity" 
condition an d possibly not even the "negativity" property (J2J). In any case the properties of the 
potential required for the validity of the various results reported below will be specified in each case. 

In the process of deriving the results presented below we also uncovered some new neat limits (such 
as those reported in the Abstract) which are as well applicable in the S-wave case and are different from 
those given in £j} . These results therefore extend those presented in Q . 

From the upper and lower limits for N# one can obtain upper and lower limits for the total number, 



(6) 



of bound states possessed by the potential V(r); the upper limit, L, of the sum in the right-hand side of 
this formula, (JHJ, is of course the largest value of I for which the potential V(r) possesses bound states. It 
is well known that the conditions © and Q are sufficient to guarantee that both L and N are finite. New 
upper and lower limits on the values of the maximal angular momentum quantum number L for which 
bound states do exist are also exhibited below, as well as new upper and lower limits on the total number 
of bound states N. Note that we are assuming, see (jfjj), to work in the (ordinary) three-dimensional 
world, with spherically symmetrical potentials. 

As in we begin below with a terse review of known results, and we then exhibit our new upper 
and lower limits and briefly outline their main features. A more detailed discussion of the properties of 
these new limits, including tests for various potentials of their cogency (compared with that of previously 
known limits), are then presented in Section m The proofs of our results are given in Section IlIII and 
some final remarks in Section [IVI 



A. Units and preliminaries 



We use the standard quantum-mechanical units such that 2m = K = 1, where m is the mass of the 
particle bound by the central potential V(r). This entails that the potential V(r) has the dimension of 
an inverse square length, hence the following two quantities are dimensionless: 



S 



dr 



-y(-)(r)' 



1/2 



(7) 



a = — max s r 

7T < ?" < oo 



-V(-\r) 



1/2 



(8) 



Here, and always below, >(r) denotes the potential that obtains from V(r) by setting to zero its 
positive part, 

y(-)(r) = V(r)9[-V(r)} . (9) 
Here, and always below, 6{x) is the standard step function, 6{x) = 1 if x > 0, 9{x) = if x < 0. 
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These quantities, S and <r, play an important role in the following. The motivation for inserting the 
(2/7r) prefactor in these definitions is to make neater some of the formulas given below. 
As for the £-wave radial Schrodinger equation, in these standard units it reads 



.2 , ^ , ^+!) 



K + V(r) 



u^K-r). (10) 



The eigenvalue problem based on this ordinary differential equation (ODE) characterizes the (moduli of 
the) £-wave bound-state energies, k 2 — ft? n , via the requirement that the corresponding cigenfunctions, 
ui(n£^ n ;r) vanish at the origin, 

ue(Ke, n ;0) = 0, (11) 
and be normalizable, hence vanish at infinity, 

lim [ut{Ki n ;r)\ = 0. (12) 

r — >oo 

It is well known that the conditions J2J and (J3J on the potential V(r) are sufficient to guarantee that 
the (singular) Sturm-Liouville problem characterized by the ODE I|1U|) with the boundary conditions 
||TJ and l|12[) have a finite (possibly vanishing) number of discrete eigenvalues ft\ n - To count them one 
notes that for sufficiently large (for definitcness, positive) values of k the solution ui(n;r) of the radial 
Schrodinger equation i|l(J|) with the boundary condition 111|) (which characterizes the solution uniquely up 
to a multiplicative constant) has no zeros in the interval < r < oo and diverges as r — > oo (proportionally 
to exp(Kr)), because for sufficiently large values of k the quantity in the square bracket on the right-hand 
side of the radial Schrodinger equation IjlOl) is positive for all values of r, hence the solution ui(k; r) of this 
second-order ODE, (fTU|) . is everywhere convex. Let us then imagine to decrease gradually the value of 
the positive constant n so that the quantity in the square bracket in the right-hand side of IjlUII becomes 
negative in some region(s) (for this to happen the potential V{r) must be itself negative in some region(s), 
this being of course a necessary condition for the existence of bound states), entailing that the solution 
U£(k; r) becomes concave in that region(s). For some value, say k = Kii, the solution ui(ni t i; r) may then 
have a zero at r = oo, namely vanish as r — > oo (proportionally to exp(— Kt t i r)), thereby satisfying the 
boundary condition at infinity <|12|) hence qualifying as a bound-state wave function and thereby entailing 
that Kg -l is the (modulus of the) binding energy of the first (the most bound) £-wave state associated with 
the potential V(r). If one decreases k below Kt t i, the zero will then enter (from the right) the interval 
< r < oo, occurring, say, at r = r^ifa) (namely ug [k; n t \ (k)] — with < r^i(«;) < oo), since the 
effect of decreasing k, by decreasing the value of the quantity in the square bracket in the right-hand side 
of the radial Schrodinger equation (|ll)fl . is to make the solution ug{n;r) more concave, hence to move 
its zeros to smaller values of r (towards the left on the positive real line < r < oo). Continuing the 
process of decreasing k, for k = a second zero of the solution ue{n; r) may appear at r = oo, entailing 
that this solution, u^K^^r), satisfies again the boundary condition at infinity 1)12(1 . hence qualifies as a 
bound-state wave function, implying that k\ 2 is the (modulus of the) binding energy of the next most 
bound £-wave state associated with the potential V(r). The process can then be continued, yielding a 
sequence of decreasing (in modulus) binding energies n\ n with n — 1,2, ...,N^. Correspondingly, the 
solution ug(n]r) of the radial Schrodinger equation IjlUI) characterized by the boundary condition (|llfl 
shall have, for K£ !rl _i > K > K£,m n — 1 zeros in the interval < r < oo. The process of decreasing the 
parameter k we just described ends when this parameter reaches the value k = 0, and it clearly entails 
that the number of zeros ri i7l (0) of the zero-energy solution u^(0; r) of the radial Schrodinger equation (pTj|l 
characterized by the boundary condition (|11|) coincides with the number Ni of bound states possessed 
by the potential V(r) (namely: ui [0; r£ jTl (0)] = with < r^i < 77.2 < . . . < r^jv f -i < ^.iVj < 00). 

Hence in the following - as indeed in jj] - in order to obtain upper and lower limits on the number 
N# of £-wave bound states we focus on obtaining upper and lower limits on the number Ni of zeros of 
the zero-energy solution ui(0;r) of the radial Schrodinger equation (|10[) characterized by the boundary 
condition (| 1 II) (for notational simplicity, these zeros will be hereafter denoted as z n , and the zero-energy 
solution of the radial Schrodinger equation 1101) characterized by the boundary condition (|llfl as u(r), 
namely u(r) = ue(0;r) with u{z n ) = 0, < z% < ... < zjsr e < 00: see Section 1111)1 . Let us moreover 



emphasize that here, and throughout this paper, we ignore the marginal possibility that the potential 
V(r) under consideration possess a "zero-energy" bound state, namely that u(r) vanish as r — > oo (for 
£ > 0), or tend to a constant value in the S-wave case; namely, we assume zn t < oo, because keeping this 
possibility into account would force us to go several times into cumbersome details, the effort to do so 
being clearly out of proportion with the additional clarification gained. 

In the following subsections we briefly review the known expressions, in terms of a given central potential 
V(r), of upper and lower limits on the number N& of ^-wave bound states, and also on the maximum 
value, L, of the angular momentum quantum number for which bound states do exist, as well as on the 
total number N of bound states, see iJBJ; and we also present our new upper and lower limits on these 
quantities. 

Before listing these upper and lower limits let us note that an immediate hunch on the accuracy of these 
limits for strong potentials may be obtained via the introduction of a (dimensionless, positive) "coupling 
constant" g by setting 



V(r) = g 2 v(r), 



(13) 



where v(r) is assumed to be independent of g, and by recalling that, at large g, Ni grows proportionally 
tog 0, 



Ni ~ g as g — ► oo, 



indeed Q 



dr 



-y(-)(r) 



1/2 



as g — > oo. 



(14a) 



(14b) 



Here, and always below, we denote with the symbols ss respectively ~ asymptotic equality respectively 
proportionality. 

The analogous asymptotic behaviors of L and of N read 



L 



9, 



(15a) 



indeed ]2c 



L max < r 

0<r<oc 



and 



indeed ll£ 



N 



3tt 



-y(->(r) 



dr r 



1/2 



as g — * oo, 



N ~ q as q — > oo 



-y(-)(r) 



3/2 



as <7 — > oo. 



(15b) 
(16a) 
(16b) 



Limits defined in terms of global properties of the potential (i. e. 

the potential) 



involving integrals over 



In this subsection we only consider results which can be formulated in terms of integrals over the 
potential V{r), possibly raised to a power, see below. We firstly review tersely known upper and lower 
limits on the number Ng of bound states, as well as known upper and lower limits on the maximum value 
L for which bound states do exist; and we then provide new upper and lower limits for both Ni and L. 

The earliest upper limit of this kind on the number Ni of £-wave bound states is due to V. Bargmann 
(and then also discussed by J. Schwinger 22|), and we hereafter refer to it as the BS£ upper limit: 



BS£: N e < 



1 



2£+l 



dr i 



-y(-)(r) 



(17) 



5 

Remark. In writing this upper limit we have used the strict inequality sign; we will follow this rule 
in all the analogous formulas we write hereafter. Let us repeat that in this manner we systematically 
ignore the possibility that a potential possess exactly the number of bound states given by the (upper 
or lower) limit expression being displayed (which in such a case would have to yield an integer value), 
since this would correspond to the occurrence of a "zero-energy bound state" (in the S-wave case) or 
a "zero-energy resonance" (in the higher-wave case) - a marginal possibility we believe can be ignored 
without significant loss of generality. 

Since the right-hand side of this inequality, iJTTJ), grows proportionally to g 2 (see $l'3\i ) rather than g 
(see l(T3|ll as g diverges, for strong potentials possessing many bound states this upper limit, l(T7|l . is 
generally very far from the exact value. It clearly implies the following upper limit on L: 



1 1 



BSL: L<L { +1 = -- + - / ,/,-,- 



OC 



V ( -\r)} . (18) 



The right-hand side of this inequality also grows proportionally to g 2 (see rather than to g (see 

(JT3J)) as g diverges, hence for strong potentials possessing many bound states this upper limit, ((TBj) . is 
also generally very far from the exact value. The limit BS£ is however best possible, namely there is a 
potential V(r), 

N e 

V(r)=gR- 1 Y^a n 5(r-/3 n R) (19) 

71 = 1 

(with an appropriate assignment of the 2Ni dimensionless constants a n and (3 n , depending of course on 
g, on t and on Ne) that possesses Ne £-w&ve bound states and for which the right-hand side of (|17f) takes 
a value arbitrarily close to Ne- 

The next upper limit we report is due to K. Chadan, A. Martin and J. Stubbe and we denote it 
as CMS. It holds only for potentials that satisfy the monotonicity condition and it reads (see 0): 




CMS: Ne < S+l- yi+ {-J £(£+1). (20) 

A less stringent but neater version of this upper limit, which we denote as CMSn, reads 

2/ 4- 1 

CMSn: N e < S + 1 — . (21) 

Clearly this inequality, (|21f) . entails the following neat upper limit on L, which we denote as CMSL: 



CMSL: /. • /. ( x:sl . \. (22) 



A more stringent but less neat upper limit on L, which we do not write, can of course be obtained from 
the CMS upper limit lf2"U|) . 

The next upper limit we report is immediately implied by a result due to A. Martin [20|, and we denote 
it as M£. It reads 



M£: N e < 



I POO 

drr 2 V e ( ±r) drV$(r) 



1/4 



(23) 



with V^ e g(r) being the negative part of the "effective £-w&ve potential" (see (tTUl ) 



(24) 
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Finally, the last two upper limits of this type we report are due to V. Glaser, H. Grosse, A. Martin 
and W. Thirring and again to K. Chadan, A. Martin and J. Stubbe an d we denote them as 
GGMT and CMS2. The first of these upper limits reads 



GGMT: N e <{2£+l) 1 - 2p C p I — 



with 



: V { -\r) 



(25a) 



Cn 



p p T 2 (p) 



(25b) 



and the restriction p > 1. This upper limit GGMT is however always characterized by an unsatisfactory 
dependence on g as g — > oo (see 't : the right-hand side of (|25a|) is proportional to g 2p with p > 1 
rather than to g (see <|13ll and l|14|0. hence it always yields a result far from the exact value for strong 
potentials possessing many bound states. 
The second of these upper limits reads 



CMS2: N e < {2£+l) 1 - 2p C p / — 

J o r 



with 



C p =p(l- P ) 



p-i 



r 2 V{r)\ 



(26a) 



(26b) 



with the restriction 1/2 < p < 1, and it is valid provided the potential is nowhere positive, see @, and 
moreover satisfies for all values of r, < r < oo, the relation 



[r^ 2p [-V(r)] 1 - p }' <0. 



(27) 



Note that the right-hand side of l)26al) features the correct power growth proportional to g, see (|13|) and 
H14f> . only if p = 1/2, in which case the condition l|27l) is equivalent to the condition but the upper 
limit CMS2 then reads (see and J7J) 



N e < 2- 3 / 2 ttS, 

hence it is analogous, but less stringent (since 2~ 3 / 2 7r = 1.11 > 1), than the limit CC, 

CC: N e < S, 



(28) 
(29) 



(see eq. (1.4) of of course this limit, valid for S-waves, is a fortiori valid for all partial waves - 
albeit clearly not very good for large £, and moreover off by a factor 2 in the asymptotic limit of strong 
potentials, see (|13|l and l|14b|) - indeed the main motivation for, and achievement of, the research reported 
in was just to provide upper and lower limits to Nq that do not have this last defect). 

Let us turn now to known lower limits, always restricting our consideration here to results which can 
be formulated in terms of integrals over the potential V(r), possibly raised to a power, see below. 

Only one result of this kind seems to be previously known 0, QJ , and we will denote it as C£. It states 
(without requiring any additional conditions on the potential other than @ and (@J) - and even these 
conditions are sufficient but non necessary) that 



C£: Ne > 1- - max 

2 7T 0<a<oo 



dr 



2C 



-aV(r)(l) 



-2f 



(30) 
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In this formula, and hereafter, the notation min [x, y] signifies x if x < y, y if y < x. Let us now assume 
that the equation 



a 



V(r) (-) (31) 



admits one and only one solution, say r = R(a) (and note that validity for all values of r of the "mono- 
tonicity condition" 10) is sufficient to guarantee that this is indeed the case) , so that the lower limit (|30[) 
can be rewritten as follows: 

11 f f R ^ ! fr\ 21 f°° f r \-2i s 

'R(a) 



N e > — + - max I dra- x (-\ - draV(r)(-) I, (32) 

2 7ro<a<oo W \ a J J R(a) \aJ 



where of course r = R(a) is the solution of (|3"T)l . It is then easy, using (|3"T|l . to calculate the maximum 
in the right-hand side of this inequality (note that the first integral inside the braces in the right-hand 
side of the above inequality, l|32|) . is elementary) and to obtain thereby the following lower limit, that we 
denote here as Cfti: 

with the radius p defined to be the solution of the following equation: 

p V(p) =(2£ + l)J dr (^) 21 V(r). (33b) 

This lower limit C£n presents the correct dependence on g, see l|13[) and (I14|l . since clearly p does not 
depend on g. It is best possible, and the potential that saturates it has the form [(J 



V(r) = -g 2 BT 2 [—) for Q<r<aR 1 (34a) 



V(r) =0 for r>aR, (34b) 
with R an arbitrary (of course positive) radius and a a dimensionless constant given by 

+ OS* < 1/2. (35) 

Let us now present our new upper and lower limits on the number of bound states possessed by the 
central potential V(r). All these results are proven in Section ITTll 

We begin with a new upper limit on the number No of S-wave bound states, possessed by a central 
potential V(r) that features the following properties: it has two zeros, V(r±) = (with r_ < r + ), it is 
positive for r smaller than r_, negative for r in the interval from r_ to r+, again positive for r larger 
than r+, 

V(r) > for 0<r<r_, (36a) 
V{r) < for r_ < r < r + , (36b) 



V(r) > for r + < r < oo. 



(36c) 



Note that we do not exclude the possibility that the potential diverge (but then to positive infinity) at the 
origin; so, for the validity of the result we now report, the condition JIJ need not hold, and indeed even 
the condition J2J can be forsaken, provided the potential does vanish at infinity, V(oo) = 0. Indeed one 
option we shall exploit below is to replace the potential V(r) with Vi !C s(r), see (I24[) . and to thereby include 
in the present framework the treatment of the £-wave case. On the other hand the assumption that the 
potential have only two zeros and no more is made here for simplicity; the extension to potentials having 
more than two zeros is straightforward, but the corresponding results lack the neatness that justifies their 
explicit presentation here (we trust any potential user of our results who needs to apply them to the more 
general case of a potential with more than two zeros will be able to obtain easily the relevant formulas 
by extending the treatment of Section ITTT|) . Let us also note that the following results remain valid (but 
may become trivial) if r_ = or r + = oo. 

We denote as NUL1 ("New Upper Limit no. 1") this result: 



NULL N <l + -Ur+ 



dr 



-V<->(r) 



1/2 



(37) 



It can actually be shown (see subsection IIII Gil that this upper limit NUL1 is generally less cogent than 
the upper limit NUL2, see Q44|l below; but it has the advantage over NUL2 of being simpler, and for this 
reason it is nevertheless worthwhile to report it separately here. 

Let us now report a new lower limit on the number Nq of S-wave bound states that holds for potentials 
that satisfy the same conditions (|3fcil> . and that we denote as NLL1 ("New Lower Limit no. 1"). It is 
actually a variation of the lower limit CI, see (|3"(T1) . and it reads: 



NLL1: N > -1 + - max 

7T 0<a<oo 



dr 



(o-Saf-^-^r)])} 



(38) 



Of course, a less cogent but perhaps simpler version of this lower bound reads 



NLLln: N > -1 



dr 



-V ( 



(r)}}{ 



max 



-V 



(-) 



(r)]} 



-1/2 



(39) 



(it clearly obtains from NLL1 by setting a = {max [— V*~ )(r)l } "'")■ 

These new upper and lower limits become relevant to the number Ng of £-wave bound states possessed 
by the central potential V(r) via the replacement in the above inequalities, H37|l and i|39[) . of V(r) with 
Ve te g(r), see l|24|l . Note that, for a large class of central potentials V(r) satisfying the conditions and 
(@J, this effective £-wave potential V£ i0 fi{r), especially for £ > 0, is indeed likely to satisfy the conditions 
(see 



F (r) = V(r) + 



e(£+l) 



> for < r < 



,.(0 



(40a) 



Vt )0ff (r) = V(r) + 



1(1+1) 



< for 



< r < r 



(40b) 



VW»0 = V( r ) + ^ > for r + ] <r<oo, 



(40c) 



required for the validity of the upper and lower limits NUL1 and NLL1, see (|37|) and (|38() . We denote 
the new upper and lower limits obtained in this manner as NUL1£ and NLLlnf: 



NULW: Ni < 1 



(r 




dr |V(r)| 



1/2 



(41) 



NLLlm? 



Ne > -1 + - 

7T 



.(«) 
+ 



dr |V(r)| 



max 



V(r) 



£(£ + !) 



-1/2 



(42) 



Next, we report new upper and lower limits on the number Nq of S-wave bound states somewhat 
analogous to those given in but applicable to nonmonotonic potentials. As above, we restrict for 
simplicity our consideration to potentials that satisfy the conditions lj36|l . We do moreover, again for 
simplicity, require the potential V(r) to possess only one minimum, at r — r m - m : 



V(r) > for < r < 



(43a) 



V(r) < 0, V'(r) < for r_ < r < r min , 



(43b) 



V{r) < 0, V'(r) > for r min < r < r+, 



(43c) 



V(r) > for r+ < r < oo. 



(43d) 



We denote these new upper, respectively lower, limits on the number Nq of S-wave bound states as NUL2, 
respectively NLL2: 



NUL2: N < l + | + -L log 



M 



(44) 



. . . . 3 •S' 1 , 

NLL2: A^ > 1 log 

2 2 2tt 8 



M 



where is of course defined by J7J and 

M = miii 



-y(-)(p),-yv->( g ) 



with the two radii p and g defined as the solutions of the following equations: 



dr 



-y(-)(r) 



1/2 



(45) 



(46) 



(47) 



-y(-)( r ) 



1/2 



7T 

2' 



(48) 



and with the additional condition (which might rule out the applicability of these limits to potentials 
possessing very few bound states, but which is certainly satisfied by potentials that are sufficiently strong 
to possess several bound states) 



P < r„ 



< 



(49) 



As already mentioned above and explained in subsection IIII Gl the upper limit NUL2, see (|44|l . is 
generally more cogent than the upper limit NUL1, see (13711 . but it requires the additional computation 
of the two radii p and q. 
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Again, as above, new limits (hereafter denoted NUL2i? respectively NLL2£) on the number Ne of ^-wave 
bound states possessed by the central potential V(r) are entailed by these results via the replacement of 
V(r) with V^eff(r), see l)24[l. so that the relevant formulas read as follows: 



Vi,es(r)>0 for 0<r<r^\ 



(50a) 



^eff(r) < 0, VU(r) < for r<? < r < r% a , 



V e<eS (r) < 0, VL s {r) > for r^ n < r < rf, 



V«, e ff(r) > for r^ } <r < oo; 



(50b) 
(50c) 
(50d) 



i r+ 

NUL2^: N e < 1 + - / dr 

7T 



1/2 1 

+ 27 



log 



^£.cft'V min) 



M 



(51) 



3 i r+ 

NLL2^: N e > — + - / dr 

2 7T .L(f) 



1/2 1 

-2^ l0g 



V £,cS I' min J 



M 



(52) 



1/2 



7T 

2' 



(53) 



+ 



dr 



1/2 



r — mm — ^ ' 



7T 

2 ; 



(54) 
(55) 



M = min 



(56) 



Let us now report another new lower limit on the number Ne of ^-wave bound states applicable to 
nonmonotonic potentials. As above, we restrict for simplicity our consideration to potentials that satisfy 
the conditions 1)36(1 . We do moreover, again for simplicity, require the potential V(r) to possess only one 
minimum, at r — r m j n , see (|43|) . We denote it by the acronym NLL3s: 



NLL3s: N e > -1 + - I dr 



- y w j '^ iog \ y(-)(^(-)( S ) )~^ iog y ' (57) 

where the radius p is defined by l|47l) and s is an arbitrary radius (of course larger than p, s > p). The 
value of s that yields the most stringent limit is a (or the) solution of the equation 



sV'(s) = 4s \V(s)\ 3/2 +4£V(s) 



(58) 
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(since clearly for this value of s the potential V(r) is negative, V(s) — — |V(s)|, in this formula V(s) 
could be replaced by V^'(s) with the condition r_ < s < r+, see 

A neater, if marginally less stringent, version of this lower limit NLL3s, which we denote as NLL3, 
reads as follows: 



NLL-'S: N t > v- - log I - 

7T \p 



(59a) 



where 



-S 



4tt 



log- 



[V ( -Hrmin)] 2 ' 

V(-)(p)V(-)(q) 



(59b) 



Here of course p is defined as above, see l|47fl . while q is defined by the formula l|48ll . and of course S is 
defined by 10. Here and below, see l|59(l and l|60|l. as well as in all subsequent formulas involving both 
p and q, see (|4"7I) and l(lS). we always assume validity of the inequality q > p, as is indeed generally the 
case for any potential possessing enough bound states. [These results are also valid if the potential has 
only one zero or no zero at all, and even if the derivative of the potential never vanishes; in this latter 
case r m i n in (|5Tjl and i|59b|l must be replaced by p] . 

Clearly this lower limit, NLL3, see l|59l) . implies the following new lower limit L^ L3L on the largest 
value L of the angular momentum quantum number I for which the potential V(r) possesses bound states 



(entailing of course that for 
bound state): 



< L 



(-) 

NLL3L 



the potential V(r) does certainly possess at least one £-wave 



NLL3L: L > L 



(-) 

NLL3L 



of course with p, q and v defined by (|47H , 14811 and J59 
integer part. 




(60) 



Here of course the double braces denote the 



C. Limits defined in terms of local properties of the potential (not involving integrals over the 

potential) 

In this subsection we report a new lower limit on the number of £-wave bound states Nt, which depends 
on the potential only via the quantity a, see (jHJ). Note that, perhaps with a slight abuse of language, we 
consider (see the title of this section) the quantity a to depend only on local properties of the potential, 
since to calculate it only the value(s) of r at which the function 2V(r)+r V'(r) vanishes must be identified. 
We also provide, in terms of this quantity a, new upper and lower limits on the largest value L of I for 
which the potential V(r) possesses bound states. 

The lower limit on Nt, which we denote NLL4, holds provided the potential V(r) satisfies, for all 
values of r, < r < oo, the inequality J5J, which as we already noted above is automatically satisfied by 
monotonically nondecreasing potentials, see It takes the neat form 

NLL4: N e > — + -. (61) 

This lower limit features the correct power growth, see ()14ajl . as g (see l|13jl) diverges, and it is best 
possible, being saturated by the potential (|34|l with 135|l . The analogy of this lower limit NLL4, see l|61|) . 

with the lower limit C£n, see (13311 . is remarkable; note that, since obviously a > {2/t:)p\V{p)\ 1 ^ 2 (see 
©), this new limit, NLL4, would always be more stringent than C^n, were it not for the additional factor 
1/2 multiplying a in the right-hand side of the inequality l|61l) (in comparison to l|33al) ). 
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This result, l|61l) . clearly entails the following new lower limit L^ll 4L on the largest value L of I for 

which the potential V(r) possesses bound states (entailing of course that for £ < £n~ll4L the potential 
V(r) does certainly possess at least one €-wave bound state): 

NLL4L: L > i« L4L = ||i(<7 - 1)}} . (62) 

Note that this lower limit features as well the correct power growth, see Q150 . as g (see C21) diverges, 
and is best possible, being saturated by the potential (|34() with i|35|) . 

Let us recall that a somewhat analogous upper limit L^g on the largest value L of £ for which the 

potential V(r) possesses bound states (entailing of course that for £ > lQ' the potential V(r) certainly 
does not possess any £-wave bound state), which we denote as ULL, reads 



ULL: L<L$> = U^(w<r-i)U. (63) 

[Indeed, it is an immediate consequence - via a standard comparison argument, see below - of the 
well-known fact that the solution u(r) characterized by the boundary condition u(0) = of the ODE 
r 2 u"{r) + cu(r) = features a zero in < r < oo only if the real constant c exceeds 4, c > 4j. 

D. Limits defined in terms of comparison potentials 

The results reported in this subsection are directly based on the elementary remark that, if VW(r) < 
V^ 2 \r) for all values of r, < r < oo, then the number of ^-wave bound states associated with the 
potential V ^ (r) cannot exceed the number of ^-wave bound states associated with the potential 

VW{r),NP>Nf?\ 

Let V(r) satisfy the negativity condition and let H^\r) be the "additional" (^-dependent) potential 
defined as follows (see Section llTlj) : 

2 



1 



with A an arbitrary nonnegative constant, A > 0. There holds then the following limits on the number 
Nt of ^-wave bound states possessed by the potential V(r): 

N e >{{XS}} if H ( x e) (r)>0 for < r < oo, (65) 



N e <{{XS}} if H^>{r)<0 for < r < oo, (66) 

where of course S is defined by Q and the double braces denote the integer part. Note however that, 
for higher partial waves (£ > 0), the lower limit, is applicable only to potentials that vanish at the 
origin (r = 0) at least proportionally to r M and asymptotically (r — > oo) no faster than r ~ 4 ( e + 1 )- while for 
S- waves (£ = 0), the upper limit is only applicable to potentials that vanish asymptotically proportionally 
to r~ p with (see ©) 2 < p < 4. 

Note in particular that (the special case with £ — and A = 1/2 of) this result implies that, for any 
potential V(r) that satisfies, in addition to the negativity condition J3J, the inequality 



there holds the following new lower bound on the number No of S-wave bound states: 

No > U I 

As can be easily verified, this lower limit is for instance applicable to the (class of) potential(s) 

,2 , „ ^ a -2 



V(r) 



9_ 

R 2 



(5)' 



exp 



/r \/3 
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(68) 



(69) 



where a and f3 are two arbitrary positive constants, a > 0, (3 > 0, that satisfy the following condition: 

a(3>0 2 + 1. (70) 

It then yields the explicit lower limit 



No > 



A 2 * r / a 

7T/3 



2/3 



(71) 



In particular, when a = 2 and (3 = 1, we obtain the lower limit No > {{2<7/-7r}} on the number of S-wave 
bound states for the exponential potential which simplifies and improves the lower bound found in our 
previous work (see eq. (2.13) of Ref. |j|). 

The particular case of the special potential V(r) that yields via l|64|) H^\r) = is investigated 
elsewhere |j . 



E. Limits of second kind, denned in terms of recursive formulas 

In this section we exhibit new upper and lower limits on the number of bound states, defined in terms 
of recursive formulas which are particularly convenient for numerical computation. We call these limits 
"of second kind," following the terminology introduced in Ref. It is possible, following Q, to derive 
such limits directly for the number Ng of €-wave bound states possessed by a central potential V(r) 
having some monotonicity properties, via a treatment based on the ODE (|191fl (see Section llllfl and 
utilizing the potential l|34|) that, as discussed in Section UTTI trivializes the solution of this ODE (just as 
the square-well potential employed in Ref. £| to obtain this kind of results trivializes the S-wave version 
of this ODE). But the results we obtained in this manner, including the precise monotonicity conditions 
on the potential V(r) required for their validity (although, as a matter of fact, the simple monotonicity 
condition Q would be more than enough for the validity of the upper limit), are not sufficiently neat, nor 
are they expected to be sufficiently stringent, to warrant our reporting them here. We rather focus on 
the derivation, using essentially the same technique employed in Ref. [4J, of new upper and lower limits 
on the number No of S-wave bound states possessed by a (nonmonotonic) potential V(r) that has two 
zeros, V(r±) = 0, and that is positive for r smaller than r_, negative for r in the interval from r_ to r + 
with only one minimum, say at r = r m j n , in this interval, and is again positive for r larger than r + , see 
(I43|l . Indeed, as already noted in subsection II Bl one can then obtain new upper and lower limits on the 
number Ne of ^-wave bound states by replacing the potential V(r) with the effective potential Ve, e ff(r), 
see Q24JI. since such a potential, for a fairly large class of potentials V(r), does indeed satisfy the shape 
conditions mentioned above: see 1|5U|). 

To get the new upper limit one introduces the following two recursions: 



(up.incr) 



^(up,incr) 



7T 

2 



(up,incr' 



) 



-1/2 



ith 



(up,incr) 
r 



(72) 



(up, deer) 



(up, deer) 



V(r 



(up, deer) 
3 



-1/2 



with 



(up, deer) 



(73) 



that define the increasing respectively decreasing sequences of radii r 



(up.incr) 



respectively 



(up, deer) 
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both 



starting from the value r m i n at which the potential attains its minimum value, see <|43(l . Now let j( u P' lncr ) 



be the first value of j such that r 



(up,i 



exceeds or equals r 



(up.incr) 

J(up,i„cr) . 



+ • 

, ~ (up.incr) 

<r+< r^J) , 



and likewise let j( u P< dccr ) 



be the first value of j such that r 



(up, deer) 
3 



(up,decr) (up.decr) 

'j(up,docr) Zi ' J(up,dccr)_ r 



(74) 

becomes smaller than, or equal to, 

(75) 



The new upper limit of the second kind (ULSK) is then provided by the neat formula 

1 r 



ULSK 



N < 



j(up,i 



j(up. 



deer) 



l 



/^(up.decr) \ 

y ./(up.dccr) J 



(76) 



Here 9{x) is the standard step function, 6{x) = 1 if x > 0, 9{x) — if x < 0, and of course r^p d d c ™] is the 
"last" (smallest) radius yielded by the recursion i)73|) , see (I75|) . [Of course the "(9-term" in the right-hand 
side of this formula, (|76|l . is not very significant, at least for potentials possessing many bound states, 
namely just when this upper limit is more likely to be quite cogent, see Section HI] . 

To obtain a lower limit one must instead define the following increasing respectively decreasing se- 

r i.. (lo,incr) ,. n (lo.dccr) 

quences ot radii respectively : 



(lo,incr) 



(lo,incr^ 



7T 

+ 2 



V(r 



(lo,incr) \ 



-1/2 



with 



(lo,incr) 
r 



> r- 



(77) 



(lo,dccr) 



n (lo,dc 



7T 

2 



V(r 



(lo,decr) \ 



-1/2 



with 



r, 



(lo,decr) 




< r+. 



Now let j( lo . incr ) be the first value of j such that - (lo ' incr) 



r 



(lo,incr) 



ncr)_! 



3 



exceeds or equals r„ 



< r y 

min _^ ' 



(78) 



(79) 



and likewise let j( u P. dccr ) b e the first value of j such that r ( uP:dccr ) becomes smaller than, or equal to, 



(lo,dccr) . 

r j(io.d C cr) S r n 



< r 



(lo,dccr) 



J(lo,docr)_ 1 - 

The new lower limit of the second kind (LLSK) is then provided by the neat formula 



LLSK 



No > \ (j (l °' ir 



r) + j(lo,dccr) _ H \ lf 



where the parameter H vanishes, H = 0, provided cither 



V 



(lo,incr) 

'jOo.incr) 



does not exceed r 



(lo,decr) 

j(lo.d C cr). 



(lo,i 



and r 



(lo,incr) 



i i i (lo.dccr) 

does not exceed r\ ( x 0idecx) , 



J) < r 

^(lo,incr) 



(lo,dccr) 

./(lo.docrJ.JJ 



^1 - r j(i°.d<=, 



/ (lo.incr) 
y J(to,incr)_ 1 



< 



mcr ) 



V r 



> 



^flOjdccr) 

(lo,decr) 



(80) 

(81) 
and 



V 



7(lo,docr)_ 1 

1, otherwise 



' J(lo,incp)_ 1 U" 1 " DAVCCU ' j(i . d ccr) J ' j(lo.incr) _ 1 ^ 1 J(lo.dccr) J ^ IS Unity, H 

[Anyway this term does not make a very significant contribution, at least for potentials possessing many 
bound states, when this lower limit is more likely to be quite cogent, see Section [H]. Note moreover 



that, in the recursions l|77|l respectively (|78|) . the starting points, 



(lo,incr) 




respectively 



restricted by inequalities; of course interesting results will obtain only by assigning 



but not exceedingly, close to r_, and r t 



(lo,decr) 




(lo,dccr) 

"o 

(lo,incr) 

"o 



are only 
relatively, 



relatively, but not exceedingly, close to r + [to get some 



understanding of which choices of these parameters, 



(lo,incr) 



and 



(lo,dccr) 



, are likely to produce more 



cogent results, the interested reader is referred to the proof of the lower limit given in Section IIIII 



of course numerically one can make a search for the values of these parameters, r 



(lo,i 



respectively 



(lo,dccr) 



' 



that maximize the right-hand side of l|81[). starting from values close to r_ respectively 
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F. Limits on the total number of bound states 



Clearly if N* ' respectively provide lower, respectively upper, limits on the number Ng of £-wave 
bound states, and likewise L^~\ respectively L' + ', provide lower, respectively upper, limits on the largest 
value L of the angular momentum quantum number I for which the potential V(r) does possess bound 
states, it is plain that the quantities 

= N fiW) , (82a) 

where 



N{L) = J2(2£+1) N^>, (82b) 

1=0 

provide lower respectively upper limits, 

JV(-) < N < N {+ \ (83) 

to the total number N, see of bound states possessed by the potential V(r). Hence several such limits 
can be easily obtained from the results reported above. 

Remark. There is however a significant loss of accuracy in using these formulas, (|82|l and (|83[1 . to obtain 
upper or lower limits on the total number of bound states N. Note that of course the upper limit, L, of the 
sum in the right-hand side of i|82bfl must be an integer, but after the sum has been performed to calculate 
N(L), see <|82b[) . it gets generally replaced by a noninteger number, respectively L^~\ to evaluate the 
upper respectively lower limit A^ + ' respectively N^~\ see (|82al) and (|83[1 . Let us illustrate this effect by 
a fictitious numerical example. Suppose we were able to prove, say, the lower limit Nt > =jp — I. We would 
then know that N > ^, Ni > ^, N 2 > |, N 3 > §, iV 4 > * entailing N > 5, Ni > 4, N 2 > 3, N 3 > 2, 
N 4 > 1 hence we could conclude that there are at least 55 bound states (5 + 12 + 15 + 14 + 9 = 55), N > 55. 
But via the above procedure we would infer that — ^ and N (L) = | (L + 1) (26 + 21 L - 4L 2 ) 

entailing > 37.23 hence N > 38. This is a much less stringent (lower) limit. Clearly, due to 

the round off errors, a lot of information got lost. This defect can be remedied, but only marginally, 
by inserting in the expression N(L) the best value of yielded by the above fictitious lower bound, 
namely = 4, since we then obtain > 38.3 hence N > 39. And the analogous calculation via 

(18211 and (|83|l from an hypothetical upper limit Ni < ^ — £, which clearly entails iVo < 4, N% < 3, 
N2 < 2, A^3 < 1 and hence N < 30, yields again a less stringent result, namely the upper limit N < 37 if 
L<+) = ±2 i s use d. This limit can be slightly improved, namely N < 35, if the integer part of L^ + ' =3 
is used. In the following we have tried to take care of this problem - to the extent possible compatibly 
with the goal to obtain simple explicit formulas. 

In the next section we illustrate the remark just made by computing firstly, via the upper and lower 
limits NUL2^ and NLL2£, two sets of integers and 7V £ (+) such that < N e < N ( e + \ and by then 

evaluating upper and lower limits on the total number TV of bound states, see (|83|l . via the standard 
formula H82b|) with L replaced, as it were, by 00, the sum being automatically stopped by the vanishing 
of the summand. The upper and lower limits obtained with this procedure will be called NUL2N and 
NLL2N respectively. 

Anyway in this subsection some results obtained via l|82|l and 183|) are reported, namely those we 
believe deserve to be displayed thanks to their neat character. But firstly let us tersely review the upper 
limits on the total number of bound states N previously known (we did not find any lower limits on N 
in the literature). 

A classical result, the validity of which is not restricted to central potentials, is known in the literature 
as the Birman-Schwinger upper bound 0iH3|, and we denote it as BiS. It reads as follows: 



BiS: N <7 ±^ [*** ^-\™^™ , (84) 



(4tt) 5 



\n - r 2 \ 
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implying, for central potentials, 

P OO P oo 

BiScentral: N<- drmV^(ri) / dr 2 r 2 V^ ] (r 2 ) log 
2 Jo Jo 

This upper limit, however, is proportional to g 4 (see l|130 ) rather than g 3 (see Ijlfijl ). hence it provides a 
limit much larger than the exact result for strong potentials possessing many bound states. 

A simple upper limit, that we denote BSN, can be obtained from the BS£ upper limit, see l|17f) : it reads 



n + r 2 



n - r 2 



(85) 



BSN: N < {{I}} 



(86) 



with 



dr r 



-y(")(r) 



(87) 



This upper limit is also proportional to g A rather than g 3 . 

An upper limit that does not have this defect and that is also valid for potentials that need not be 
central was obtained by E. Lieb [l8j- We denote it as L: 



L: N < 0.116 / d 3 f 



-y(-)(f) 



3/2 



(for the origin of the numerical coefficient in the right-hand side of this formula, we refer to the original 
paper [Ig)- For central potentials it reads as follows: 



Lcentral: N < 1.458 



dr j 



-y<-)(r) 



3/2 



(89) 



(the numerical coefficient in this formula is of course obtained by multiplying that in the preceding 
formula by 4-7r; for other results of this kind, none of which seems however to be more stringent than 
those reported here, see ||). 

Let us end this listing of previously known results by reporting the upper limit on the total number of 
bound states N obtained [ll| by inserting (|21[1 and (|22H in 182|l . As entailed by its origin, it only holds 
for monotonically nondecreasing potentials, see and (J2J). We denote it as CMSN: 



CMSN: N < 



12 



S 3 + 3S 2 



3- 



2tt 



(90) 



with S defined by Q. 

We did not obtain any new upper limit on the total number ./V of bound states sufficiently neat to be 
worth reporting. We report instead a rather trivial upper limit on TV obtained via (|82[) with L replaced 

by its upper limit L$ (see JHSJ) and with N e < N and N bounded above by NUL2, see (jSJ. This 
upper limit on the total number N of bound states is therefore applicable to potentials that satisfy the 
condition (|36() . and we denote it as NUL2Nn: 



NUL2Nn: N < -(na + if \ 2 + S + - log 



M 



(91) 



with M defined by a defined by (©, S defined by and V^~\r m - m ) the minimal value of (the 

negative part of) the potential, see Ij36|l . For a monotonia potential, see |T]l. this upper limit takes the 
simpler form 



NUL2Nm: N < -{ira + l) 2 \l + S 



2n 



log 



v(p) 

V(q) 



(92) 
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where p and q are defined by i|47|) and (|48|l respectively (this result is of course obtained using, instead 
of NUL2, the analogous result valid for monotonic potentials [4|). It is remarkable that, in spite of the 
drastic approximation Ng < No used to get these two limits, they turn out, in all the tests performed in 
Section m to be more stringent than all previously known results. 

We now report two new lower limits, which recommend themselves because of their neatness, although, 
for the reason outlined above, one cannot expect them to be very stringent. 

A new lower limit, that we denote as NLLN3, on the total number N of bound states for a potential 
V(r) that satisfies the conditions l|36|) . follows from the lower limits NLL3, see (|59|) . and NLL3L, see (|60|) . 
A simple calculation yields 



with v defined by (|59b|) and 



NLLN3: N > — (2v + \)(v + A), 
6X Z 



A=ilog(? 



(93) 



(94) 



with q and p defined by l|47|) and (|48|l (we assume of course q > p, hence A > 0). 

Another new lower limit, that we denote as NLLN4, on the total number N of bound states for a 
potential V(r) that satisfies the condition (JjJJ is implied, via l)82[l. from the lower limits NLL4, see l|61|l . 
and NLL4L, see (1621) . A simple calculation yields 



NLLN4: N > 



1 



1 



(95) 



with a defined by ©■ Here, as usual, the double braces denote the integer part. This lower limit has the 
merit of being rather neat, but it grows proportionally to g 2 (see JTJJJ) rather than g 3 (see l)16fl). hence 
it cannot be expected to be cogent for strong potentials possessing many bound states. 

II. TESTS 

In this section we test the efficiency of the new limits reported in Section |U by comparing them for 
some representative potentials with the exact results and with the results obtained via previously known 
limits. For these tests we use three different potentials: the Morse potential 21] (hereafter referred to as 
M) 



(T 

a 

\R 



exp 



M: V{r) = ~g 2 iT 2 |2exp 
the exponential potential (hereafter referred to as E) 

E: V(r) =-g 2 R- 2 exp (™) 
and the Yukawa potential (hereafter referred to as Y) 

Y: V(r) = -g 2 {rR)- 1 



(5-)]}' 



ex P( -- 



(96) 



(97) 



(98) 



In all these equations, and below, R is an arbitrary (of course positive) given radius, and g, as well as a 
in H96|l . are arbitrary dimcnsionless positive constants. 



Tests of the limits on the number of bound states Nt 



The first potential we use to test the new limits is the M potential l|96|l. This is a nonmonotonic 
potential for which the number Nq of bound states for vanishing angular momentum is known; we indeed 
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consider for this potential only the I = case. [We do not test the GGMT and CMS2 limits with this M 
potential since, from their incorrect behavior when the strength g of the potential diverges, we already 
know that these limits give poor results. But, in spite of this incorrect behavior, these limits could be 
useful when there are few bound states; they are therefore tested below with the E and Y potentials, in 
cases with £ > 0] . 

The exact formula for the number of S-wave bound states for the M potential is 



N 



(99) 



Note that it is independent of the value of the constant a. 

For this potential, the limits NUL2 and NLL2, see (|44fl and (|45p. can be computed (almost completely) 
analytically: 



NUL2 



NLL2 



No 



No 



< 9 
> 9 



^-logs + 1, 

Z7T 

1 1 3 

— log s , 

2vr 6 2' 



with s = min(2y — y 2 , 2x — x 2 ) and x, y solutions of 



7r - yj2y-y 2 - 2 arcsin f |j 
3in (I) 



2 arcsin i 



7T 

2? 

7T 

2^ 



(100) 
(101) 

(102) 
(103) 



The calculation of the cutoff radii p and q, see (|47|) and (|48(l . cannot be evaluated analytically. But one 
can compute upper and lower limits, p < p and q > q, on these radii by using only the attractive part of 
the potential in the definition l|47|) and l|48l) of p and q. When p and q are used in place of p and q we 
obtain the (marginally less stringent) limits (denotes as NUL2s and NLL2s) 



NUL2s : 



< £f+— log 

Z7T 



NLL2s : N > g 



Z7T 



4(z 2 - 1) 
4(z 2 - 1) 



1, 

3 
2' 



(104) 
(105) 



with z = 8g/ir. As mentioned in Section [I] validity of the inequalities p < r m j n < q is required in order 
to use the NUL2s and NLL2s limits; this leads to the restriction g > 7rv / 2/(8(v / 2 — 1)) = 1.34. 
The NLLl limit l|38|l takes for the M potential the simple form 



NLLl: N > 0.672 g - 1. 



(106) 



The other new limits cannot be tested with this potential: the upper limit NUL1 and the limits of the 
second kind are not applicable because r + = oo, the lower limit NLL4 is only applicable to monotone 
potentials, and the lower limit NLL3 coincides with NLL2 for I = 0. 

The previously known limits applicable to this potential (note that the CMS upper limit is only 
applicable to monotone potentials) take the form: 



BS^ 



A r o<2.g 2 (a + log2+-); 



(107) 



Ml 



No < V2g(a 2 + (3 - 21og2)a + 1.901) 



1/4. 



(108) 



C£: No > -i + lgT(a), 

Z 7T 



T(a) 



max 

0<7<1 



exp(a) 



- exp(2a) - a/1-7 2 + y log _ { 



1 + VI 



(109a) 
(109b) 
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FIG. 1: Comparison between the exact value l|99|l of iVo (ladder), the upper limits BS<? l|l()7fl (long dash), 
ifTHHl (spar e dot ) and NUL2s (short dash) and the lower limits 0,1 fTDTty (dash-dot), NLL2s ifTH^jl (solid) 

and NLL1 itTUrJl (dash-dot-dot) for the M potential |@ (all with I = and a = 1). 



To obtain this limit we set a — R/(-fg) in (|30fl with £ = 0. Note that for a > uq = 1.386, with 
4exp(ao) = exp(2ao), this lower limit C£ is trivial because T(a) is then negative). Note that, for 
a = log 2, T(a) reach is maximal value: T(a) — 1.055. The factor multiplying g in the right-hand side 
of l|lU9afl coincides which that multiplying g in the right-hand side of the NLL1 lower limit (|106[) : so for 
this particular value of a, this C£ limit is slightly more stringent (of course for all values of g) than the 
NLL1 lower limit. For all other values of a, there exists a value g a such as for all g > g a , the NLL1 limit 
is more stringent than the C£ limit. For example, for a — 1, the NLL1 limit yield more cogent results 
than the C£ limit as soon as g > 5.22, namely as soon as the number of bound states is greater than five. 

Fig. ^displays these limits as a function of g. Note that the limits depend on a while the exact result 
does not. We tested the results for the a = 1 case (not a = 0, in oder to have a nonmonotonic potential). 
It is clear from this figure that the generalizations to nonmonotonic potentials of the results obtained in 
Ref. , namely the limits NUL2 and NLL2, are quite cogent. This remains true even for large values of 
g : for instance, when the exact number Nq of bound states is 5000, these upper and lower limits restrict 
its value to the rather small interval [4996, 5003]. In this case the BS£ upper limit exceeds 1.5 10 8 , the 
M£ upper limit only informs us that N < 10307, the lower limit C£ that N > 2879 and the lower limit 
NLL1 that N > 3359. 

The second test is performed with the E potential (|97|l . The exact number Ni of bound states for this 
potential is computed by integrating numerically (|191fl with 77(0) = and 77(00) = N^tt, see Section ITTT1 
The upper limit NUL1£ reads 

NULW: N t < 1 + ~V%+ - x- jV [exp(-a;_) - exp(-x+)] - l[i + 1) - -^j j , (110) 
where x± are the two solutions of 

£(£+ 1) = g 2 x 2 ± exp(-x±). (Ill) 
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The NLLln£ limit reads 



NULln^ : N t > -1 + ~ |.g 2 [exp(-z_) - exp(-z+)] - £(£ + 1) (J- - -^j J \V^§\-^ 2 , (112) 

where is the minimal value of the effective potential (|24|l . The NUL2^ and NLL2^ limits can be 

written as follows: 



NUL2£: N t < -F(g,£; ar_, x+) + ^- log 

7T 27T 



NLL2£: iV £ > £; a;_, x+) - — log 

7T 27T 



T/min 



A/ 



T/min 



M 



where 



and 



F(g,£;a,b) = / — a/ g 2 x 2 exp 



(-as) -*(*+!), 



M = min[|^, eff (p)|,|^, eff (g)|] 



where p and g are solutions of 



F{g,£-x.,^) = ^, F(g,£;^x + ) = \- 



The lower limits NLL3 and NLL4 take much simpler forms: 



NLL3: N t >-g- ^log(-) -§--log 

7T Z7T V 7T / 2 7T 



log a; 



l0g(l - 3C) 



with x = 7r/(4g), and 



NLL4 : TV* > 



2.9 



1 



7re(2£+l) 2' 



The previously known limits are found to be: 



BS^ : N e < 



9 



21 + r 



(113) 
(114) 

(115) 
(116a) 
(116b) 

(117) 
(118) 
(119) 



CMS : N t < — + 1-\1 + 



A£{£+1) 



(120) 



M£ : N e < (AB) 1/4 

where A and B are given by 

A = g 2 [exp(-x-)(x 2 _ + 2x- + 2) - exp(-x + )(x+ + 2x+ + 2) 

1 1 
x - x + , 



B = g 2 [exp(-ar-) - exp(-x+)] - + 1) 



(121a) 

e(£+l)(x + -x-), (121b) 
(121c) 
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TABLE I: Comparison for the E potential 1971 between the exact number Nt of bound states, various upper and 
lower limits on Nt previously known and new upper and lower limits on Nt, for several values of g and I. 



LLSK NLL3 NLLln^ NLL2^ Ex NUL2£ BS^ CMS M£ GGMT NXJLll ULSK 
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3 


3 


3 
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21 
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21 


12 
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1 
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9 
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(-; 
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1 3 

-LO 
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c; 
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4 


4 


■-> 

o 




« 
o 




Id 




'HI 

O -L 


-LO 


q 
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2 
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2 


3 


4 


13 


13 


7 


7 


7 


5 


IS 


3 


7 


6 


6 


7 


9 


10 


46 


21 


18 


43 


25 


11 




9 


2 
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2 


4 


4 


17 


17 


9 


8 


9 


5 


24 


4 


10 


8 


8 


10 


12 


13 


64 


28 


24 


60 


33 


14 




12 


4 
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4 


5 


6 


23 


23 


13 


11 


12 


7 


29 


5 


13 


9 


9 


13 


15 


16 


76 


34 


29 


71 


40 


17 




15 


4 





5 


4 


6 


7 


27 


28 


15 


13 


13 


8 


35 


6 


16 


11 


11 


16 


18 


19 


94 


41 


35 


88 


49 


20 




18 


6 





6 


6 


7 


8 


33 


33 


18 


16 


16 


9 


40 


7 


18 


12 


13 


18 


20 


21 


106 


47 


40 


100 


55 


22 




21 


6 





6 


6 


8 


9 


37 


38 


20 


18 


18 


10 



GGMT: N e < g 2 ?(2l + l) (1 ~ 2p) C ^ 2p \ 



(122) 



with C„ defined by (|2 



CMS2 : Ne < g 2 P(2l + l)^"^ 



(123) 



with C p defined by (|2 



C£ : N f > 



2.9 



tt(2£ + 1) 



yexp 



( D 2' 



(124) 



where y is the solution of yexp(— y) = (2£ + l)y 2e T(l — 21, y). 

Comparisons between the various limits and the exact results are presented in Tabled The BS£ limit 
gives poor results when g becomes large but becomes slightly better as i grows. The CMS gives better 
restrictions when I is small but behaves like the BS^ limit when t grows. The M^ limit overestimates the 
number of bound states by a factor 2 when I is small; it is no better for larger I, yet better than the BS^ 
and CMS limits. The GGMT limit (with, in each case, the optimized value of the parameter p, see (112211 ) 
gives similar results to those yielded by the BS^ limit when I is small and becomes better and equivalent 
to the M£ limit for larger values of i. The results obtained with the CMS2 limit are uninteresting hence 
not reported: indeed, the values of p which minimize the value of the limit are either p = 1/2 for small 
values of I (in which case this limit is analogous but less stringent than the CC limit, see H29I) 1 ). or p = 1 
for larger values of £ (and this yields the BS^ limit). The new limits NUL2^ and NLL2i? clearly yield the 
most stringent results. The NLLln£ lower limit only yields cogent results for large values of the angular 
momentum. The NLL3 lower limit works reasonably well for small values of £ but becomes poor for 
higher values of the angular momentum. The limits of the second kind ULSK and LLSK yield similar 
results to those given by the NUL2^ and NLL2^ limits. Note that the arbitrary radii 7-^ lo,mcl - ) and r ( lo,dccr ' 
have been chosen to optimize the restriction on the number of £-wave bound states. Finally, the results 
obtained with the C-f'n and the NLL4 lower limits are not reported because they are very poor. These 
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limits give Ne > 1 for small value of I and Nt > for large value of I. This defect comes from the 
presence of the factor l/(2t + 1) which for instance implies that this lower bound becomes three times 
smaller when I go from to 1 while the actual number of bound states Nt decreases generally only by 
one or two units. 

The last test is performed with the Y potential 1)98(1 . The exact number Nt of bound states is again 
computed by integrating numerically (|191|l with 77(0) = and 77(00) = Nfir, see Section ITTT1 
The NUL1£ limit takes the form 

1/2 

g 2 £ + dx e ^~ x) - l{£ + 1) (J- - J- U , (125) 

where x± are the two solutions of the following equation 

1(1 + 1) = g 2 x± exp(-x±). (126) 

The NLLln£ limit reads 

NLLln£ : N £ > -1 + 1 jg 2 £ dx ^ v{ ~ x) - 1(1 + 1) (J- - -L) j |^|" x / 2 , (127) 

where V^™g is the minimum value of the effective potential 124|) . The NUL2£ and NLL2^ limits can be 
written as follows: 



1 

-( 

n 

1 

-( 

7T 



1 


log 


T/imn 


2^ 


M 


1 


log 


T/min 


2^ 




M 



3 
2' 



(128) 
(129) 



where 



and 



/ — y/g 2 xexp( 

Ja X 



G(g,i;a,b)= — ^/g 2 xey^(-x) - £(£ + 1). (130) 



M = min [\V e<eS (p)\, \Vt >efi (q)\} (131a) 

where p and q are solutions of 

G( ff ,£;x_,|) = |, G(g,l;±x + ) = 1 (131b) 
The lower limits NLL3 and NLL4 take somewhat simpler forms: 

NLL3; „ (> ^ 9 _^^)_|_(I±^), 

where y and x are defined by erf(y) = \J n/8/g, erf(x) = 1 — -J n/8/g, and 

NLL4: N <> l^fc+T) -* (133) 

The previously known limits take the form: 

BS^ : N t < (134) 
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TABLE II: Comparison for the Y potential 198H between the exact number Nt of bound states, various upper and 
lower limits on Nt previously known and new upper and lower limits on Nt, for several values of g and I. 



LLSK NLL3 NLLlnf NLL2^ Ex NUL2£ BS^ CMS M£ GGMT NUL1I ULSK 
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53 


40 
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39 


46 
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20 
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41 
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23 
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47 
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10 


45 


54 


20 


20 


22 
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2 / 4£(£+l) 

CMS: N t <2J~g+l-dl+ (135) 

: N e < (AB) 1/4 (136a) 

where A and B are given by 

A = g 2 {exY>(-x-)(l + X-)-exp(-x + )(l + x + ))-l{£+l){x + -x_), (136b) 

B = g 2 f + dx eM ~ x) -M + l)(- -); (136c) 

J x X \X- X+J 

GGMT: N e < g 2p {2£ + \)^p) ■ ( 137 ) 

CMS2: N t <g 2p {2t+l) {1 - 2 ^ 9l£M ■ ( 138 ) 



Ct: ^>-^L_Vv«p(-|)-i, d39) 

where y is the solution of exp(— y) = (21 + l)y 2£ T(— 21, y). 

Comparisons between the various limits and the exact results are presented in Table [H] The char- 
acteristics of the various limits for the Y potential are analogous to those commented above for the E 
potential. Here again the NUL2£ and the NLL2£ limits are the most effective ones, being indeed fairly 
stringent for all the values of g and I considered, and the limits of the second kind ULSK and LLSK also 
give quite stringent limits. 
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B. Tests of the limits for the value of L 



In this subsection we test various limits on the largest value L of the angular momentum quantum 
number I for which the potentials E and Y do possess bound states. In this article we only obtained new 
lower limits on L. Indeed neat upper limits on L cannot be extracted from the new upper limits on Ng 
presented in Section [I] But we will now see that the "naive" upper limit L^g is quite good indeed 

better than the previously known upper limits and £cm SL , see C^l 1 an d JUJl- 

The first test is performed with the nonsingular E potential (|97(l . The lower limit NLL3L gives 



L 



(-) 

NLL3L 



with 



2g _l_ 

7T 2lT 



log 



1 -x 



A = -log 

7T 



logx 



log(l - x) 



and x — 7r/(4g). The lower limit NLL4L gives 



3 
2' 



(140a) 

(140b) 
(140c) 



L 



NLL4L 



1 ( 4g 



2 \-rre 



- 1 



(141) 



The previously known upper limits on L can also be obtained analytically: 

1 , 



r(+) _ 
BSL 



(142) 



L 



(+) _ 

CMSL — 



- 2 {49-D 



(143) 



-(+) 



1 (4g 



- 1 



(144) 



A comparison between these limits and the exact results (computed numerically, as indicated in the 
preceding subsection) is presented in Fig. Except for the BSL upper limit, we have the correct linear 

behavior in g as discussed above, the best result being clearly provided by the naive upper bound . 



-(+) 



is of 



Indeed this result appears hardly improvable, because the error introduced by this upper limit L eS 
at most one unit (at least for this example, as well as the following one, see below). 

The second test is performed with the singular Y potential (again the exact result can only be computed 
numerically). The lower limit NLL3L gives 



with 



L 



(-) 

NLL3L 



{{- 



(145a) 



x 



A = 



2, 

-log 

7T 



;t 2n 
x 
V 



y 2 i , 



(145b) 
(145c) 
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FIG. 2: Comparison between the exact value of L (diamond), the upper limits BSL 11421 (solid), CMSL (11431 (long 
dash) and fllH (short dash) and the lower limits NLL3L flltl (dash-dot) and NLL4L IT3T1 (dash-dot-dot) 
for the E potential (1571 . 



where y and x are defined by erf(y) = a, erf(a^) = 1 — a, a = yj (tt/8)/ g. The lower limit NLL4L gives 

-(-) ff 1 , 9 



L 



NLL4L 



2 ?rVe. 

The previously known limits on L can also be obtained analytically: 

(+) 



L 



BSL 



(146) 



(147) 



-(+) 

^CMSL 



(148) 



(149) 



The comparison between these limits and the exact results is presented in Fig. [3J Here again, the best 
result is obtained from the naive upper bound L^g for which the error is, as indicated above, of at most 
one unit. 



C. Tests of the limits for the total number of bound states TV 



This subsection is devoted to test the limits on the total number of bound states N. We will not test 
the BiS limit due to its bad behavior at large g. We do however test the BSN limit which yields the same 
incorrect behavior but is simpler to compute. 



2G 
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FIG. 3: Comparison between the exact value of L (diamond), the upper limits BSL 11471 (solid), CMSL 11481 (long 
dash) and L { c p ( IHSty (short dash) and the lower limits NLL3L (THESl (dash-dot) and NLL4L iTHroT (dash-dot-dot) 
for the Y potential 19811 . 



The first test is performed with the E potential (|97|l . Again, the exact result can only be calculated 
numerically. 

The new upper limit NUL2Nm, see l|92|l . takes the simple form 



8\e / V 7T 7T 7T 



The new lower limit NLLN3 reads 



N 



>^(2 (4ll3l) 2 + 74ll3L + 6), 



(150) 



(151) 



with ^jjll3L an< ^ v Si yen by equations (|140a|) and l|140b|) ; and the new lower limit NLLN4 is given by 



NLL3L 

(|9"5j) with a = 4g/(e7r 

The other, previously known limits read as follows: 



BSN : N < ~g 2 (g 2 + 1), 



(152) 



Lcentral : N < 0.864 g 3 , 



(153) 



CMSN 



N < 1.698 (g 3 + 2.3562 5 2 + 1.116 g + 0.1473). 



(154) 



The BSN and the CMSN limits can be improved: instead of using the limit on L provided by these 

limits (iggL an d -^cmsl)> we can use * ne best u PP er limits L^t' H63|l: the BSN and CMSN limits obtained 
in this manner are called here improved BSN and CMSN limits: 



improved BSN: N < —g 



(155) 
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FIG. 4: Comparison between the exact value of N (diamond), the upper limits BSN 11521 (solid), CMSN 11541 
(long dash), Lcentral (ffixfl l (short dash), improved BSN iffiaBTl (dot), improved CMSN iTHpflt (sparse dot), NUL2Nm 
iTHafjl (dash dot dot) and the lower limits NLLN3 itTCTl (dash dot) and NLLN4 (ladder) for the E potential 
1971 . The black and the white circles correspond respectively to the NUL2N and NLL2N limits. 



improved CMSN: N < 0.5202 (g 3 + 2.179 g 2 + 1.726 g + 0.4806) (156) 

Fig. ^presents a comparison between the various limits and the exact result. It shows that the limits 
on the total number of bound states which can be expressed in a neat form are not very stringent. 
[Indeed, the best results are yielded by the upper limit NUL2Nn which is obtained using only a limit on 

the number of S-wave bound states, No, and the simple limit on the maximal value L of I for which 
bound states do exist]. There are at least three reasons for this. First, most of the limits do not contain 
the appropriate functional of the potential (as identified by the asymptotic behavior at large g of N, see 
(|16bp 'l: only the Lieb limit Lcentral, see l|89|l . features the correct form, but the numerical factor is not 
optimal indeed too large (by approximately a factor 7). The second reason is that for every value of I, 
there is a round-off error introduced by the limit; to obtain the limit on the total number of bound states 
we sum all these errors. The third reason is that to be able to make the summation over the values of i 
we must have an explicit dependence of the limits on £, and this entails that we cannot use some of the 
limits we found; in particular we cannot use the new upper and lower limits NUL2 and the NLL2, which 
are quite stringent, to obtain a neat formula. But we can use them and compute upper and lower limits 
on Nf, then sum all these contributions to obtain upper and lower limits on TV, the sum being stopped 
when N« + ' is smaller than 1 and N,~' is negative (see subsection II F|) . We call NUL2N respectively 
NLL2N the upper respectively lower limits on the total number of bound states N obtained (from NUL2 
respectively NLL2) via this (inelegant) procedure. Fig. 0] shows, for 4 values of g, that these limits are 
quite stringent. 

The second test is performed with the Y potential. Again, the exact total number of bound states is 
computed numerically, as indicated above. 

The new upper limit NUL2Nm, see l|92(l reads 

NUL2Nm: N<~(^ + l) \2g\[^ + — ^ + - log - + 1 ) , (157) 

8 \V e / \V7r 7T 7r y I 
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where y and x are defined by erf(y) = y/Tt/8/g, erf(ir) = 1 — \/ir/8/g- The new limits NLLN3 and 
NLLN4 are: 



NLLN3: N > "-(2 (4ii 3L ) * + + 6 J , (158) 

with 4ll3L an< ^ v Si ven ^ equations (|145fl and (|145bjl : while the lower limit NLLN4 is given by (|95|l 
with a = 2g/{y/eir). 

The other previously known limits read: 

BSN : N < l -g 2 {g 2 + 1), (159) 



Lcentral: iV< 0.703.g 3 , (160) 



CMSN: 7V< 3.3422( ff 3 + 1.88.g 2 + 0.711 5 + 0.075). (161) 



1 2 ^2 5 



improved BSN: N < -g 2 \ + 1 ) (162) 



improved CMSN: N < 0.4924 (g 3 + 2.237 g 2 + 1.781 g + 0.5078) (163) 

A comparison between the various limits and the exact numerical results is presented in Fig. El in 
analogy to the case of the E potential and with analogous conclusions, see above. 

III. PROOFS 

In this section we prove the new results reported in Section [i] Because we tried and presented those 
results in Section [I] in a user- friendly order, the proofs given below do not follow the same order, due to 
the need here to follow a more logical sequence. To provide some guidance we divided this section into 
several subsections, but we must forewarn the reader that a sequential reading is essential to understand 
what goes on. 

A. Tools 

The starting point of our treatment is the well-known fact (see Section HJ that the number N# of (£- 
wave) bound states possessed by the potential V(r) coincides with the number of zeros, in the interval 
< r < oo, of the function u(r) uniquely defined (up to an irrelevant multiplicative constant) as the 
solution of the zero-energy ^-wave radial Schrodinger equation 



u"(r) 

with the boundary condition 



u(r) (164a) 



u(0) = 0. 



(164b) 
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FIG. 5: Comparison between the exact value of N (diamond), the upper limits BSN 11591 (solid), CMSN 11611 
(long dash), Lcentral 11601 (short dash), improved BSN 11621 (dot) and improved CMSN 11631 (sparse dot), 
NUL2Nm IT571 (dash dot dot) and the lower limits NLLN3 IT55I (dash dot) and NLLN4 (ladder) for the Y 
potential 1981 . The black and the white circles correspond respectively to the NUL2N and NLL2N limits. 



To get an efficient handle on the task of counting these zeros (or rather, of providing upper and lower 
limits on their number) it is convenient to introduce a new dependent variable rj(r) related to u(r) as 
follows: 

l-U(r)} 1/2 cot [„(r)] = f(r) + (165) 

Here we moreover introduce two new functions: a function f(r), which we reserve to assign at our 
convenience below with the only proviso that it be finite in the open interval < r < oo, and nonnegative 
throughout this interval, 

f(r) > for < r < oo; (166) 

and a function U(r), which might or might not coincide with the potential V(r) (see below) but that 
(unless we explicitly state otherwise) we require to be finite in the open interval < r < oo , to satisfy 
(at least) the properties (see ©, © and (@J) 

-U(r) = \U(r)\ , (167) 
lim \r 2+£ U(r)} = 0, (168) 

r — >oo 



lim [r 2 - £ U(r)] = 0, (169) 



and to be related to the potential V(r) as follows: 



U(r) + W(r) = V{r) + ^ ± l) , (170) 
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where we still reserve the privilege to assign the function W(r) at our convenience (a possibility will be 
to set U(r) = V(r) hence W(r) — £(£ + l)/r 2 ; but it shall not be the only one, see below). Of course the 
function r](r) (as well as u(r)) depends on £, although for notational simplicity we omit to indicate this 
explicitly, and this remark may as well apply to the other functions, f(r), U(r), W(r), introduced here 
and utilized below. 

It is then easy to see that the function r\ (r) is uniquely characterized by the first-order nonlinear ODE 
(implied by ljTF5|) with lfT7Ujl . and I|l{j4a|) l 

n '{r) = \U(r)\ 1/2 + \U(r)\- 1/2 {lf(r)] 2 -f(r)-W(r)} sin 2 [„(r)] 

-{[^Pir)]}- 1 U'(r) + f(r)}^n[2 V (r)}, (171) 

with the boundary condition (implied by 11651) with J164b|) and l|169|) ) 

r,(0) = 0. (172) 

We of course assume the function r/(r) to be continuous, disposing thereby of the mod(7r) ambiguity 
entailed by the definition 1165(1 . 

It is now easy to sec that this function rj(r) provides a convenient tool to evaluate the number of zeros 
Ni of u(r). Indeed, if we denote with z n the zeros of u(r), u(z n ) — 0, ordered so that 



= z < zi < . . . < z Nl 



(173) 



(see Section nj, since clearly (see l|171|) l whenever ij(r) is an integer multiple of tt the derivative rf{r) is 
nonnegative, rj'(r) > 0, we may conclude (see H165|l ) that 



ry(z„)=n7r, n = 0, 1, . . . , N e , 



with 



(n — l)ir < rj{r) < nn for z n _i < r < z n , n 
It is moreover plain that, provided (see pfi5[l ) 

\U(r)\ 1/2 



lim 

r — ^oo 



fir) 



u' (r) 
u(r) 



o, 



there holds the asymptotic relation 



77(00) = NiTT, 



(174) 
(175) 

(176) 

(177) 



and that this asymptotic value is approached from above. [To prove the last statement one sets, in 
the asymptotic large r region, rj{r) = Ngn + e(r) with |e(r)| <C 1, and uses the asymptotic estimates 
sin 2 (A^7r + e) w e 2 , sin [2(Nen + e)] ~ 2e to rewrite the ODE (|1 71|> in the asymptotic region as follows 
(recall that U(r) is nonpositive, see l|167)) 1: 



e'(r) « {-U(r)f 2 + |^ - 2 /(r)} e(r). 



(178) 



One can then replace the approximate equality sign rj in this formula with the equality sign = and 
integrate the resulting linear ODE, obtaining 



e(r) = [[/(r)! 1 / 2 / dx exp 



dyf(y) 



(179) 



where R is a (finite) integration constant, and this formula (valid in the asymptotic, large r, region, where 
of course r > R) shows that e(r) is indeed positive]. 
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To get a more detailed information on the behavior of the function r)(r) we make the additional 
assumption that the assignments of the auxiliary functions f(r) and W(r) guarantee validity of the 
following inequality: 

[/Wf - f'(r) - W(r) > (180) 

(it would be enough for our purposes that this inequality be valid only at large values of r - but for 
simplicity we assume hereafter its validity for all values of r, < r < oo). It is then clear from the ODE 
(|171fl satisfied by rj(r) that, wherever rj(r) takes a value which is an odd integer multiple of 7r/2, namely 
at the points r = b n such that 

» 7 (6 tt ) = i(2n-l)7r ) n = 1, 2, . . . ,N i} (181) 

its derivative rj'(r) is nonnegative, T]'(b n ) > 0. Hence we may complement the information provided by 
(jl74f) with that provided by this formula, l|181|) , and moreover replace the information provided by (|175fl 
with the following more detailed information: 

(n- 1)tt < v(r) < ^(2n- 1)tt for z„_i < r < &„, n = 1, . . . (182) 

(2n - 1)| < »j(r) < riTr for h„ < r < z„, n = 1, . . . , N e , (183) 

which of course also entails that the points z n and b n are interlaced, 

z < h < Z\ < b 2 < ■ ■ ■ < Ziv f -i < bjf t < ZN e < oo. (184) 

And note in particular that these formulas entail the following important inequality (implied by the non 
existence of 6jv^+i), valid for all values of r, < r < oo: 

< v(r) < (n ( + I) 7T. (185) 



2 



It is moreover clear that the maximum value of f](r), 



fj = max [r](r)] , (186) 

0<r<oo 



is actually attained in the interval zn ( < r < oo, and that it lies in the range 

N e ir < i) < ( N t+^\ 7T. (187) 
We are now in the position to derive the new upper and lower limits reported in Section [I] 

B. Proof of the lower limits NLL4 

We begin by proving the lower limit NLL4, see H(jlf) . To this end we assign as follows the functions 
f(r) and W(r); 

f{r) = 1 -, W(r) = ^±^-, (188) 

entailing that the left-hand side of the inequality (|180fl vanishes, that 

U(r) = V(r), (189) 



that the definition (|165f) now reads 
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\V(r)\^ cot [„(r)] = ~ + ^, 



and, most importantly, that l|171|) reads 



^ ) = |F(r)|1/2 -{iM + ^ sin[2,?(r)1 



(190) 



(191) 



Here we are of course assuming the potential V(r) to satisfy the condition J2J. 

Before proceeding with the proof, let us note that, for the potential l|34l) . the second term in the 
right-hand side of this ODE, l|191fl . vanishes, hence one immediately obtains 



7](r) 



fJ 



2£+l 



2(+l 



r \2£+i 

r) 



Hence (see ty!77\i ). for the potential 



Nt = 



ga 



21+1 



(192) 



(193) 



where as usual the double braces signify that the integral part must be taken of their contents. This 
observation implies that the upper and lower limits which obtain (see below) by massaging the last term 
in the right-hand side of the ODE i|191|) are generally best possible, being saturated by the potential (|34|l 
(if need be, with an appropriate choice of the parameter a, see Section 0. 

To prove Hfilj) with (jHJ we now introduce the auxiliary function r?io(^) via the ODE 



Vi(r) = W(r)\ 1/2 [-^^ + ^} mo(r) (194) 

with the boundary condition 

mo(0) = 0. (195) 

We then assume the potential V(r) to satisfy, in addition to J3J), the condition (0. It is then plain (see 
(tT^Hl with (|T72l . (|TM|> with (|T35|) . and ©) that, for all values of r, < r < oo, 

mo{r) < v(r). (196) 

[Indeed (|194|l obtains from 1191|l via the replacement sin(x) =>■ x, and for positive x, sin(x) < x; hence 
f]\ (r) can never overtake rj(r) because, at the overtaking point, a comparison of l|194(l with ljl91|l entails 
r l\ { r ) — r l'{' r )i which negates the possibility to perform the overtaking]. But the linear ODE (|194[) with 
(I195fl can be easily integrated to yield (recalling (0)) 

hence we conclude (see (|T3oT) . (THJ7I) and ©) that 

* * 2T27TT) ' (198) 



and via l|187|l this entails the lower limit NLL4, see (|61fl . which is thereby proven. 



C. Proof of the upper and lower limits NUL2 and NLL2 
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Let us now prove the upper and lower limits NUL2 and NLL2, see l|44|l and (|45|l . on the number No of 
S-wave bound states possessed by the central potential V(r), of course under the assumption that this 
potential satisfy the conditions l|43|l . The main tool of the proof is the same function rj (r) as defined in 
the preceding subsection IIII Bl which is therefore now defined by the formula (see H190|i ) 



\V{r)\ 1 / 2 cot[r ) (r)} = ^, 
uyr) 



and satisfies the ODE (see dU) 



rf(r) = \V(r)\ 1 ^-^Lsm[2r 1 (r)}. 



(199) 



(200) 



Note that this is just the function rj(r) that provided our main analytical tool in Q ; however the conditions 
Ia|) and (|43b|) (which clearly imply V(r_) = 0) entail now, via H199|) . the condition 

(201) 

hence it has no 
(202) 



ry(r_) = 0, 

as well as the fact that u(r) is concave in the interval < r < r_ (see (|lG4a[l and 
zero in that interval, hence 

r_ < bi < z\ 



(see (|174f) . (I181f) and H199I) ). Likewise the fact that u(r) is also concave in the interval r + < r < oo (see 
(|164a|) and (|43djl ). hence it has no extremum in that interval, entails 



zn -i < b No < r + . 



(203) 



To obtain the upper limit NUL2 we now integrate the ODE (|200|l from Z\ to zm -i (and note that, 
thanks to (|5U51) and (|3U3j) . as well as g3J), we can hereafter replace, whenever convenient, |^(r)| with 
-y(")(r), see @): 



r)(zN -i) - = (N -2)ir 

1 f z «o- 

~4 



ZN -1 



dr 



dr 

V'(r) _ s 



\V(r)\ 



sin [2ry(r)] . 



1/2 



V'(r) 

dr WH shl[2ri(r)] 



(204) 



The first equality is of course entailed by (|174|) . As for the second equation, note that we conveniently 
split the integration of the second term in the right-hand side of i|200|) in two parts. The properties 
H43 bfl . (|43cl) (as well as the obvious fact that |sin(cc)| < 1), allow us to majorize the right-hand side of 
this equation, (|204|l . We thereby get 



(N - 2) 7T < 



z iV -l 



dr 



-y(-)(r) 



1/2 



dr 



V'( r ) 1 p- 1 ^ V'(r) 



(205) 



hence, 



{No - 2) 7T < 



/ dr 


-V { -\r) 


1/2 









+ i log- 



[V(~\r min )Y 



Vl-)(zi)V(-)(z No -i) ' 



(206) 



We now need to find quantities p and q, defined only in terms of the potential, such that p < Z\ and 
q > 2jv -i an d also such that \V(p)\ < \V(zi)\ and \V(q)\ < \V(zn -i)\. Let us first consider the 
"favorable" case (which obtains for a sufficiently attractive potential): z\ < r m i n and zat _i > r m ; n . In 
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this case, we integrate (|200|) from b\ to z% and since in this interval both V'{r)J \ V(r)\ and sin [2ry(r)] are 
negative, we infer 



T7(^i) - 77(61) = - < / dr -V(~\r) 



1/2 



(207) 



hence a fortiori (see (120211 s ) 



If we define p via the formula 



r 

Jo 

f 

Jo 



dr 



dr 



-y(-)(r) 



-V(->(r) 



1/2 



> 



1/2 



(208) 



(209) 



then we conclude (by comparing l|2U8|l with 112091) ) that p < Z\. Moreover, since we have supposed 
z i < r min, we have also |^(p)| < | V(^i)|. We then integrate (|200|l from zn -i to b No , and taking 
advantage of the fact that in this interval both V'{r)J \V(r)\ and sin [2?7(r)] are positive, we infer 



v(bN ) - r)(z No -i) = 77 < 



dr 



-y(-)(r) 



1/2 



(210) 



hence a fortiori (see (120311 s ) 



7T 

2< 



dr 



-y(->(r) 



1/2 



(211) 



Analogously, if we define g via the formula 



1/2 



(212) 



we conclude that q > zn _i. Moreover, since we have supposed zn -i > r m - ln , we have also |y(g)| < 
|y(zjv -i)|- Thus if these two relations, z\ < r m i n and zn„-i > r m in, hold, we obtain 



N < 



dr 



-y(-)(r) 



1/2 



47T 



log- 



[^"Hrmin)] 2 

v(-)(p)y(-)(«) 



(213) 



where we have used the equations (|209|) and l|212|) . This last relation imply the validity of the relation 

We need now to consider the cases where z\ > r m ; n or z^-i < f* m ; n . In these cases we could have, for 
example, p < z\ and \V(p)\ > |V(zi)|. 

Firstly, suppose z\ > r min which implies that ZjVo-i > r m ; n . Now we impose a first condition for the 
applicability of the limit NUL2: p < r min . This condition is always true for a potential which possesses 
enough bound states; in practice, the limit NUL2 will be applicable only when the attractive strength of 
the potential is large enough. This condition ensures that p < z\ (this is not necessarily true, with the 
definition (|2()9|) . when z\ > r m i n ). Indeed, if z\ < r m i n , we have proved it above, and if z\ > r m i n , this is 
still true since p < r m ; n . Moreover, since z\ > r m ; n , we have | V(^i)| > |F(zAr _i)| > |V(g)| > M, with 



M = mm(\V(p)\,\V(q)\) 



-y(-)(p),-y(-)( g ) 



(214) 



This implies the validity of the relation (|44fl . 
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Secondly, suppose Zjv -i < r m i n which obviously implies z\ < r m ; n . Now we impose a second condition 
for the applicability of the limit: q > r m j n . This condition is always true for a potential which possesses 
enough bound states. This condition ensure that q > zm —i (this is not necessarily true, with the definition 
()212p . when zat _i < r m i n ). Indeed, if zn ~i > ^mim we have proved it above, and if Zjy -i < r min, this 
is still true since q > r m [ n . Moreover, since 2/v -i < fmin, we have |V(«jv — i)| > l^( z i)l > |^(p)| > M, 
with M defined again by (|214|l . This implies, via the definitions Q and 10, the validity of the relation 
(|44|l . and concludes our proof of the new upper limit NUL2. 

The proof of the new lower limit NLL2, see (|45|) . is completely analogous, except that one integrates 
the ODE HP from p to q 



dr 



-V { -\r) 



1/2 



dr 



V'(r) , 1 
V(r) J 



9 dr^ 
-1 „ V(r) ■ 



and from the inequalities (actually valid for any positive radius, see (1185(1 ) 

T)(P) > 0, 



(215) 



(216) 



77(g) < [No + - 7T, 



we clearly infer 



(217) 



Hence 



N 



7T > 



dr 



1/2 



-y(-)( r ) 

hence, via the definitions l(209|l and Q212JI of p and g, 



N 



IT > 



dr 



-y(-)( r ) 



1/2 



J V"(r) 1 
V(r) 4 
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V(r) 



log- 



[^-Hrmin)] 2 

W-)(p)W-)(g) 



(218) 



(219) 



(220) 



Note that this inequality is true in any case, provided p < r m j n < q, and of course it implies (again, via 
the definitions 10 and ©, as well as ((214(1 ) the validity of the marginally less stringent lower bound ((45(1 
(we preferred to display in Section [I] the lower bound ((4511 rather than the more stringent one implied by 
((220(1 to underline its analogy with the upper bound 144(1 ). 



D. Proof of the lower limits NLL3s and NLL3 



Let us now proceed and prove (following |4|) the new lower limits NLL3s and NLL3, see ((57(1 and 
((59(1 . The proof is analogous to the proofs of the NUL2 and NLL2 limits given in the previous subsection 
except that now instead of considering the equation ((200(1 for 77 (r) we use the equation ((191|) . To obtain 
NLL3s, we integrate the ODE l(191|) from p to an arbitrary radius s > r m ; n : 



7?(s) - r/(p) 



[ dr 


-y(-)(r)" 


1/2 1 r 


Ip 







V'(r) 
dr sin [2 V (r)} 

\V(r)\ 



V'(r) 

£ 

dr - sin [2ry(r)] . 
r 



(221) 
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The right-hand side of this last equation can be minorized (since rj(s) < (Ng + 1/2)tv and rj(p) > 0, see 
0S3) to yield 



Ni + - j it > rj(s) - r)(p) > I dr -V^\r) 



1/2 1 

~ 4 



dr 



V'(r) 1 
V(r) + 4 



V'(r) 
V(r) 



From the definition of p, see (I47|) . we finally obtain 



(N e + 1)tt > / dr [-V (_) (r) 



1/2 i 



•log- 



' [^-Hrmin)] 2 

y(-)(p)v(-)(a) 



£loe 



VP/ 
(222) 



(223) 



which coincides with the lower limit NLL3s, see l|57ll . that is thereby proven. 

To get the lower limit NLL3 we proceed as above, except that we integrate from p to q, see (I47f) and 

6EJ>- 



E. Proof of the results in terms of comparison potentials (see Section II DH 



Let us now prove the relations l|65|l and Ij66|l . We assume for this purpose that the potential V(r) 
satisfy the negativity condition but we require no monotonicity condition on V(r); we do however 
require the potential V(r) to be nonsingular for < r < oo and to satisfy the conditions, see Q and Q, 
that are sufficient to guarantee that the quantity S, see Q, be finite. 

Let us now replace the potential V(r) with V(r) + H^\r), so that the radial Schrodinger equation, 
see l|164all . read now 



u"[r) 



£(£+1) 



u(r) 



and the relation j 170(1 read now 



U(r) + W{r) = V(r) + H ( x e) (r) + 



e(e + i) 



(224) 



(225) 



Let us moreover set (see H165I) ) 



U(r) = V(r), 



(226) 



so that (see 



f(r) 



V'(r) 
AV{r) ' 



(227) 



(228) 



with the additional requirement 



flf (r) + W+H [f{r) f + f {r) = f3 | V(r )| , (229) 

where j3 < 1. [Note that we imposed in l|180|l that the left-hand side of l|229|) be positive. Actually this 
restriction, which was introduced to prove that T]'(b n ) > (see I|181|l 1. was too strong for our needs. 
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Indeed, one can verify from l|171(l that we still have rf{b n ) > provided j3 < 1]. It is easily seen that this 
entails for H^' (r) the definition (|64f> via the assignment 

/3 = 1-4A 2 . (230) 
Hence these assignments imply that the definition of ?7(r), see (|165H . reads now 

[-^)] 1/2 cot[, W ] = ^ + ^4, (231) 
AV{r) u(r) 

and, most importantly, that the equation i|171l) satisfied by r/(r) becomes now simply 

V '(r) = \V{r)\ 1 ' 2 {1 -0 sin 2 [77(f)]}, (232) 

entailing 

r/(r) {1-/3 sin 2 [r/(r)]} _1 = |^(r)| 1/2 . (233) 

Both sides of this last equation are now easily integrated, the right-hand side from r = to r = oo, 
and the left-hand side, correspondingly, from r\ = (see (|172fl . which is clearly implied by l|231|) and by 
(I164b|l l to 77(00), yielding (see 10 and l|2*3Tjj l) 

7/(00) = XSn. (234) 
It is on the other hand clear that in this case as well 

N# 7T < 77(00) < (JVg + 1) 7T. (235) 

[Indeed, while in this case the relation (|176|l does not hold and therefore neither l|177fl nor (|187|l need 
be true, the relation l|174|l is still implied by the definition l|231|) . and moreover (|233|l clearly implies 
il'( z n) > 0, entailing validity of these inequalities]. Hence (sec 1234f) ') 

N^ +H ^ ] ={{\S}} , (236) 

(V-\-H ) 

where as usual the double brace denote the integer part. In this last formula the notation A 
denotes of course the number of £-wave bound states possessed by the potential V{r) + H^\r). 

But if the "additional potential" H^\r), see l|64f) . is nowhere negative, this potential V(r) + H^\r) 
cannot possess less (^-wave) bound states than the potential V(r), hence the lower limit l|65|l is proved. 
And under the same conditions, if the function H^\r) is nowhere positive, the potential V{r) + H^\r) 
has no less bound states than the potential V(r), hence the upper limit (|6l)|l is proven. 

F. Proof of the upper and lower limits NUL1 and NLL1 

Next, we prove the new upper and lower limits NUL1 and NLL1, see l|37|) and l|38|) . To prove them we 
of course assume the potential V(r) to possess the properties (|36fl. and we set I = 0. We moreover set 

(see dinni)) 

f{f) - (237) 

and 



U(r) = 



(238) 
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where a is a positive constant, a > 0, the value of which we reserve to assign at our convenience later. 
Note that in this case our assignment for U(r) does not satisfy the conditions (|168fl and (|169[) . and that 
the definition (|165(l of rj(r) now reads 



a cot [r](r)] — 



u'(r) 
u(r) 



Consistently with these assignments we also set (see 

W(r) = V(r)+a 2 , 
and the equation satisfied by r](r), see (|171f) . now reads 

7]' (r) = a cos 2 [rj(r)] — aT 1 V(r) sin 2 [r](r)] . 
We then integrate this ODE from r = r_ to r = r + : 



dr {a cos 2 [rj(r)] + a^ 1 \V(r)\ sin 2 [r]{r)}} 



(239) 

(240) 
(241) 

(242) 



Note that we used (|36b|) . 

Now the definition l|239|) of rj(r) entails that the radii b n , see (|181|l . coincide with the extrema of the 
zero-energy wave function u(r), u'(b n ) — 0. We therefore can use l|36af) to conclude that, since the zero- 
energy wave function u(r) is concave in the interval < r < r_ , the first extremum b\ must occur after 
r_, r_ < b\, hence (see p81|l ) 



< 7j(r_) < 



(243) 



Likewise, l|36c|l entails that u(r) is concave in the interval r + < r < oo, hence the last extremum, bpf , 
must occur before r+, 6jv < r+, while of course r](r) can never reach (No + 1) ir (note that in this case 
the condition l|176[l does not hold hence the more stringent condition (|185|l does not apply). Hence 



No--) n< V {r+) < (N + 1) tt, 



(244) 



where we are of course denoting as No the number of bound states possessed by the potential V(r). 
Hence we may assert that 



(N - 1) tt < 77(7 + ) - ry(r_) < (N + 1) tt. 
From the left-hand one of these two inequalities, and Q242JI . we immediately get 

(No-1) tt< / + dr {a \V(r)\}, 

J r_ 

since the replacement of cos 2 [r](r)] and sin 2 [rj(r)] by unity in the right-hand side of 
(further) majorization. Hence 



No < 1 



a (r + — r_ ) + a 1 



dr 



-V { -\r) 



and by setting 



v-l/2 



rf?" 



-y(-)(r) 



1/2 



(245) 

(246) 
entails a 

(247) 
(248) 



we get 



iV < 1 + 



(r+ - r_) 



-V ( ->(r) 



1/2 



39 



(249) 



The result NUL1, see l|37[) . is thereby proven. 

To prove the lower limit NLL1, see (|38|) . we use the second inequality of l|245|) to get from (|242|l the 
inequality: 

poo r r -, -, 

(250) 



poo 

/ dr min 


a, a 1 


-y(")(r)" 
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(Wo + 1)tt > 77(r+) -7j(r_) > 
and this coincides with the formula l|38|l . The lower limit NLL1 is thereby proven 



G. Unified derivation of the upper limits NUL1 and NUL2 

In this subsection we indicate how the derivations of the two upper limits, NUL1 and NUL2, can be 
unified. This entails that, in this context, the limit NUL2 is the optimal one. In view of the previous 
developments our treatment here is rather terse. 

The starting point of the treatment is the ODE 



r/(r) = \U(r)\^ 2 cos 2 [„(r)] - V(r) \U(r)r^ 2 sin 2 fo(r)] + ^\ sin[2ry(r)] 

4 U (r) 



(251) 



that corresponds to (|171|l with I — 0, f(r) — 0, and W(r) = V(r) — U(r), where we always assume 
U(r) to be nonpositive, see (|167[) . but otherwise we maintain the option to assign it at our convenience. 
Clearly this formula entails 



r]'(r) < max \\U{r)\ 1/2 , -V(r) \U{r)\' 1/2 



1/2" 


+ 


U'(r) 






4U{r) 



hence 









max 







r){r 2 ) - r)(n) < dr I max \\U(r)\ 1/2 ,-V{r) |[/(r)f 1/2 



U'(r) 



AU(r) 



(252) 



(253) 



It is now clear that two assignments of U (r) recommend themselves. One possibility is to assume that 
U(r) is constant, implying that the last term in the right-hand side of this inequality vanishes: this is 
indeed the choice (|238|l . and it leads to the neat upper limit NUL1, see the preceding subsection llll Fl The 
other, optimal, possibility is to equate the two arguments of the maximum functional in the right-hand 
side of this inequality, (|253|l . namely to make the assignment i|189|) . and then to proceed as in subsection 
IIII Dl arriving thereby to the upper limit NUL2 (and note that a closely analogous procedure was used 
in albeit in the simpler context of a monotonically increasing potential). 



H. Proof of the upper and lower limits of second kind (see Section II ED 

Let us now prove the results of Section lTEl beginning with the proof of the upper limit (|76|l . To this end 
let us assume first of all that r^p d '^ ) > (see (|75|l and {7SJ), and let us then introduce the piecewise 
constant comparison potential V^ + '(r), defined as follows: 

V {+) (r)=0 for 0<r<rJP;' 



(up, deer) 

doer) : 



(254a) 



y(+)(r) = V(r^ dccl) ) for r ^ decl) < r < r ^ p ' decl) with j = j( u P' dccr ', j( u P' dccr ) - 1, . . . , 0, 

(254b) 
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V(+\r)=V(rf p ' ma >) fo 



(up, iner) Nj f (up, iner) 



< r < r 



(up, iner) 



with j = 1, . . . , J 



(up,incr) 



(254c) 



F (+) (r) = for r J(U p, i „ cr ) < r. 
It is obvious by construction (if in doubt, draw a graph!) that 

V(r) > V i+) (r), 



(254d) 



(255) 



hence, if wc indicate with Nq the number of S-wave bound states possessed by the potential (r) 
(and of course with N the number of S-wave bound states possessed by the potential V(r)), clearly 



N <K 



(+) 



(256) 



It is moreover clear from the previous treatment, see in particular l|199|) and 1200|) that we now write, 
in self-evident notation, as follows, 



V (+) (r) 



1/2 



cot 



r/ + >(r) 



ttW(r) 
u(+)(r) ' 



r/+)'(r) = ^ (+) (r) 



i/2 y(+)'( r ) 



4|V(+)(r)| 



2?7 (+) (r) 



(257) 



(258) 



that the number JVq of S-wave bound states possessed by the potential V^ + \r) can be obtained in 
the usual manner via the solution 7/ + '(r), for r > r^^ dc "2), of this ODE, i|258[l . characterized by the 
boundary condition (see 1)2571) and l|254a|l entailing u(r) = u'(Q) r for < r < 



V y, /(up, doer) J 



1/2 



cot 



-,(+) 



namely 



entailing 



tan 



„(+) /„("P.decr) \ 



/„(up.decr) \ 

y J(up, deer) J 

v V J(u| 



n (up,decr) 

J(up. door) 



r (up,dccr) x 

J(up, doer) , 
1 



„(up,decr) 

j(up, deer) 



deer) 

up, deer) 



1/2 



-,(+) 



(V up 
y 



deer) \ 7T 

up, deer) ) ~0 ' 



(259) 



(260) 



(261) 



The number iV (+J of S-wave bound states possessed by the potential V(+) (r) is then characterized by the 
inequality 

N M<i v ( + )^^ + 1 _ (262) 

[Indeed, in self-evident notation, b^ +) < r^^j}^, see l|254d|) . hence r]^ f&^c+)^ — V 
(see (ESHl and and [ft^j J = - ±) tt, see (tl%Tl) 1. 



(+) / T ,< up,incr) 

^/(up, iner) 



We now introduce another solution, 7/ ++ )(r), of the ODE l|258[l in the interval 



(up, iner) 

J(up, iner) 



„(up,decr) 

J(up, deer) 



characterized by the boundary condition 



j) 



(++) fJiP.decr) \ _ ^ 

^' J(up, deer) j ~ ^' 



< r < 



(263) 
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It is then plain, from a comparison of this "initial condition" l|263|) with l|261|) . that, throughout this 
interval, r/ ++ ^(r) > ?/( + )(r), hence in particular f rjy"^""^ ^ > rf + ^ (r^^ 

we infer a fortiori 



.mcr 



j i , hence from 



A r o (+) <^ (++) (^r.T2,)+^ (264) 
But the function 7/++) (r) can be easily evaluated in closed form, since de facto it satisfies the ODE 

r,(++)'( r ) = 7<+)(r) 1/2 . (265) 

Indeed the second term in the right-hand side of the ODE (|258|l now vanishes: inside the intervals in which 
the potential V^ + \r) is constant, see (|254ll . because its derivative V^'(r) vanishes; and at the boundary 
of these intervals, where the potential V^(r) is discontinuous hence its derivative V^'(r) features a 
delta-function contribution, because the term sin [2?7 ( - ++ ^(r)] vanishes: indeed, as can be immediately 
verified, the ODE (|265[) with the initial condition (|263[) and the piece- wise potential (I254|l . entails that 
at these boundaries, say at r = r„ with 

„(+) _ „,( u P,dccr) „(+) _ JuP.dccr) _(+) _ ^(uP.decr) 

- ' J(up, decr)_ 1 ) • • • ) ' J(up, doer) — ' 1 ) 



_(+) _ r . _(+) _ Jup.incr) (+) _ (up,incr) 

' J(up, docr) +1 — 'mill, ' J(up. decr)_|_2 _ 1 ' J("P- dccr) +3 ~ '2 > 

„(+) (up.incr) (+) _ (up.incr) 



t)_l) 'j(up, dccr) + j(up. incr) + 1 — ',/(up, incr) ) (ZOOj 

there hold the relations 

?7 (++) (r r ( + ) ) = n |, n = 1, 2, j( up - dccr ) + J< up ' incr > + 1. (267) 



This last formula entails indeed sin 



0, and moreover 



hence, via (|264|l . 



V i++) (^Ti)) = (/ (up ' dccr) + ^ (up ' incr) + l) J, (268) 

JV <j ( j( u P' incr ) + j("P.dccr) + _ pggj 



consistently with the new upper limit of the second kind, see J7BJ, in the case r^^^jr) > 0. 

If instead r^Ju P ?!w) < 0, the proof is analogous, except that the formula l)254a|l is now irrelevant, the 
"initial condition" l|261|) is replaced by 

77 (+) (0) = 0, (270) 

the "initial condition" 1|263[) is replaced by 

-1/2 



V 



( (up, deer) \ 

I ' J(up,dccr) _l J 



(0) = TT v \ v , (271) 

' 2 (up, deer) v ' 

r Jtup.dccr).! 

in the definition l|266[) of the radii r„ ' the lower index in the right-hand side is always decreased by one 
unit (i. e., r[ +) = r^; d ^ ) _ 1 and so on, entailing r^2 Pi docr)+J(up , incr) = ry££ n ™J T) ), hence (|268|l reads 

*? C++) (^toci)) = (j (up ' dccr) + J (up ' incr) ) J, (272) 



and consequently 1269(1 reads 



N < - (j(»P^«) + j(up,dccr) + ^ 
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(273) 



consistently with (|76|) in the case 



(up. deer) 

J(up, deer) 



< 0. The proof of the new upper limit of the second kind 

(|76|l is thereby completed. 

The proof of the new lower limit of the second kind, see JHTJ, is analogous, and we therefore only 
outline it here. It is again based on the construction of a piecewise potential V^~'(r) that (in this case) 

maximizes the potential V(r) and for which the number Nq 'of S-wave bound states (or rather, an upper 
limit to this number) can be computed easily in closed form. We manufacture this piecewise potential 
V^~\r) according to the following prescriptions: 



' (r) = oo for < r < r, 



(lo,i 



(274a) 



t r( — ") / \ \T ( (lo,incr)\ r (lo,incr) ~ (lo.incr) 

y\ , {r)=V\r j J 1 '\ for r)_\ < r < r) 



'j-i 



with j = 1, 2, ... , j0°. incr ) _ 1, (274b) 



y(-)( r ) = 



V 



( (lo.incr) \ r, ( (lo.decr) 



for < r < r^ (l0 ,,., r 



(lo,dccr) 



(274c) 



^-'(r) = y(r] lo ' docr) ) for rg' dccr) < r < s 

with j = j( lo . de «') - 1, j(lo,decr) _ 2, . . . , 1, 0, 



t (lo,decr) 



(274d) 



V"( )(r) = max[V(r)] for r ( lo ' dccl ' < r _ 

r->r+ 



(274e) 



It is then obvious (draw graph if in doubt!) that this potential maximizes the original potential V(r), 



V(r) < V { - ] {r), 



(275) 



hence that the number Nq -*of its bound states provides a lower limit to the number iVo of bound states 
of the potential V(r), 



N >N n 



(-) 



(276) 



But it is also clear, on the basis of the analysis given above (and leaving to the alert reader the task 
to provide the details required to turn this argument into a rigorous proof), that each of the j( lo > lncr ) -)- 
j(io,decr) m t erv als in which the piecewise potential V^(r) is negative, see (j274bT) . l|273c| and pTiHji . 
can accommodate "half a bound state" (namely it yields an increase by 7r/2 of the relevant function 
77(f)), except possibly for the central interval around r m i n , which can or cannot accommodate such "half 
a bound state" depending whether the product of the square root of the modulus of the potential V^'(r) 
in that interval times the length of that interval, 



(lo,dccr) (lo,incr) 

? "j(lo,docr)_ 1 ~ r j(lo,i„cr)_ 1 



v f r ( lo < inc1 ') ^ V /V lo ' docr ) ^ 

V ^J(lo,i„o r) _ ly ) , ' ^j(lo,decr)_ 1 J 



1/2 



(277) 



does or does not amount to no less than 7r/2 (and it easy to verify that the definition of H as given 
after l|81|) entails H = in the former case, H = 1 in the latter: see (|77ll and 178fl ). This justifies the 
expression in the right-hand side of l|81|l . except for the additional term —1 appearing there, which takes 
care of the two extremal intervals of the potential V^'(r), see (|274a|l and l|274e|) . each of which can at 
most "unbound half a bound state" (i.e., cause a decrease of the relevant r](r) by tt/2). 



IV. OUTLOOK 
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With modern personal computers the calculation of the number of bound states in a given central 
potential V(r) is an easy numerical task (especially using 1|191|) with i|172|) and as we ^ as the 

numerical computation of the corresponding binding energies and eigenfunctions. It remains however 
of interest to obtain neat formulas which provide directly in terms of the potential V(r) upper and 
lower limits for these physical quantities - as indeed demonstrated by the continued attention given to 
these problems in the recent literature 0,0, H EE EH EE EE EE EE E3- The technique used in this 
paper, and in the one that preceded it go back to the 60's (see for instance 0), yet our findings 
demonstrate that it can still yield remarkably neat, and cogent, new results. We plan to explore in the 
future the applicability of this approach to establish upper and lower limits to the energies of bound 
states, as well as to obtain upper and lower limits on the number of bound states (or, more generally, of 
discrete eigenvalues) in more general contexts, including those spectral problems on the entire line that 
are relevant for the investigation of integrable nonlinear partial differential equations, a context in which 
the number of discrete eigenvalues is generally related to the number of solitons, see for instance 0. 
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